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degree  of  Doctor  of  Philosophy  in  Chemical  Oceanography 

Abstract 

Zinc  (Zn)  stable  isotopes  can  record  information  about  important  oceanographic 
processes.  This  thesis  presents  data  on  Zn  isotopes  in  anthropogenic  materials, 
hydrothermal  fluids  and  minerals,  cultured  marine  phytoplankton,  natural  plankton,  and 
seawater.  By  measuring  Zn  isotopes  in  a  diverse  array  of  marine  samples,  we  hope  to 
understand  how  Zn  isotopes  are  fractionated  in  the  oceans  and  how  Zn  isotopes  may  be 
used  as  tracers  of  marine  biogcochcmical  processes.  Common  forms  of  anthropogenic  Zn 
had  b^Zn  from  +0.08  %o  to  +0.32  %o,  a  range  similar  to  Zn  ores  and  terrigenous 
materials.  Larger  variations  were  discovered  in  hydrothermal  fluids  and  minerals,  with 
hydrothermal  fluids  ranging  in  b^Zn  from  0.02  %o  to  +0.93  %o,  and  chimney  minerals 
ranging  from  -0.09  %o  to  +1.17  %o.  Lower-temperature  vent  systems  had  higher  bw'Zn 
values,  suggesting  that  precipitation  of  isotopically  light  Zn  sulfides  drives  much  of  the 
Zn  isotope  fractionation  in  hydrothermal  systems.  In  cultured  diatoms,  a  relationship  was 
discovered  between  Zn  transport  by  either  high-affinity  or  low-affinity  uptake  pathways, 
and  the  magnitude  of  Zn  isotope  fractionation.  We  established  isotope  effects  of 
b6,,Zn  =  -0.2  %o  for  high-affinity  uptake  and  b66Zn  =  -0.8  %o  for  low-affinity  uptake.  This 
work  is  the  first  to  describe  the  molecular  basis  for  biological  fractionation  of  transition 
metals.  Biological  fractionation  of  Zn  isotopes  under  natural  conditions  was  investigated 
by  measuring  Zn  isotopes  in  plankton  collected  in  the  Peru  Upwelting  Region  and  around 
the  world.  Seawater  dissolved  Zn  isotopes  also  reflect  the  chemical  and  biological 
cycling  of  Zn.  The  bA6Zn  of  deep  seawater  in  the  North  Pacific  and  North  Atlantic  is 
about  0.5%o,  and  the  dissolved  b6ftZn  gets  lighter  in  the  upper  water  column.  This  is 
unexpected  based  our  observations  of  a  biological  preference  for  uptake  of  light  Zn 
isotopes,  and  suggests  that  Zn  transport  to  deep  waters  may  occur  by  Zn  adsorption  to 
sinking  particles  rather  than  as  primary  biological  Zn.  The  thesis,  by  presenting  data  on 
several  important  aspects  of  Zn  isotope  cycling  in  the  oceans,  lays  the  groundwork  for 
further  use  of  Zn  isotopes  as  a  marine  biogcochcmical  tracer. 


Thesis  supervisor:  Edward  A.  Boyle 
Title:  Professor  of  Ocean  Geochemistry 
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Chapter  1 
Introduction 


Zinc  (Zn)  atoms,  regardless  of  isotope,  have  thirty  protons  in  the  nucleus  and 
therefore  have  a  similar  configuration  of  electrons  orbiting  the  nucleus.  The  chemistry  of 
Zn  is  determined  mostly  by  the  interactions  of  these  protons  and  electrons  with  each  other 
and  with  charged  particles  of  other  atoms.  The  difference  between  Zn  isotopes  lies  in  the 
number  of  neutrons  in  the  nucleus.  Slight  differences  in  nuclear  mass  will  have  a  small 
effect  on  the  strength  of  the  bonds  that  the  atoms  forms.  Slight  differences  in  the  bonding 
strength  of  difference  isotopes  can  lead  to  a  preferential  accumulation  ol  heavy  or  light 
Zn  isotopes  in  certain  bonding  environments.  Zn  isotope  ratios  thereiore  contain 
information  about  the  chemical  processes  that  have  acted  upon  Zn  atoms.  The  growth  of 
phytoplankton  in  the  oceans,  chemical  adsorption  to  sinking  particles,  the  circulation  of 
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seawater  through  hydrothermal  vents,  and  other  oceanographic  processes  involve  changes 
in  the  chemical  bonding  of  Zn.  Zn  isotopes  therefore  contain  information  about  the 
chemical  history  of  samples  and  may  be  valuable  tracers  of  chemical  oceanographic 
processes. 

1.1  Historical  overview 

1.1.1  The  development  of  MC-ICP-MS 

The  development  of  multi -col lector  inductively-coupled  plasma  mass 
spectrometry  (MC-ICP-MS)  has  increased  the  ease  and  accuracy  of  isotope  measurement 
for  trace  elements.  Isotopic  analysis  by  thermal  ionization  mass  spectrometry  (TIMS)  is, 
in  comparison,  more  time-consuming,  less  effective  at  ionizing  transition  metals  such  as 
Zn,  and  subject  to  errors  based  on  the  time-dependent  fractionation  of  isotopes  during 
evaporation  from  a  filament.  The  accuracy  of  single-collector  ICP-MS,  in  which  different 
isotopes  are  measured  in  rapid  succession,  is  limited  by  plasma  instability,  rapid 
variations  in  ionization  efficiency  on  timescales  less  than  the  timescale  of  measurement. 
Multi-collector  measurements  correct  for  this  by  simultaneously  collecting  data  on  all 
masses  of  interest. 

The  development  of  MC-ICP-MS  has  led  to  an  explosion  in  the  study  of  “non- 
traditional”  stable  isotopes.  Natural  variations  have  been  discovered  in  the  stable  isotopes 
of  Cr,  Fe,  Cu,  Zn,  Mo,  Cd,  Hg,  and  many  other  elements  (Anbar  and  Rouxel,  In  revision; 
Johnson  el  a!.,  2004). 
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1.1.2  Natural  Zn  isotopes  in  environmental  and  marine  systems 

Figure  1  presents  the  natural  Zn  isotope  measurements  that  had  been  reported  at 
the  time  1  proposed  to  undertake  this  thesis.  Data  on  Zn  isotopes  in  the  ocean  included 
measurements  on  manganese  nodules  and  marine  sediment  trap  material  (Marechal  et  al., 
2000).  Manganese  nodules  had  generally  higher  Zn  isotope  ratios  in  polar  regions  than  in 
low  latitudes,  an  effect  that  was  suggested  to  correlate  with  biological  Zn  uptake  in  the 
upper  ocean  (Fig.  2).  Sediment  trap  data  showed  a  remarkable  seasonal  signal  in  the  Zn 
isotope  composition  of  sinking  particles  (Chapter  6),  Data  showing  large  variations  in  the 
6!’<’Zn  leached  from  bulk  sedimentary  carbonates  over  the  last  1 75  ka  suggested  that  Zn 
isotopes  might  have  uses  as  a  paleotracer  for  surface  productivity  (Pichal  et  al.,  2003). 

Several  notable  papers  on  Zn  isotopes  in  marine  and  environmental  systems  have 
been  published  in  the  past  few  years.  Zinc  isotopes  have  been  measured  in  ancient 
hydrothermal  deposits  (Mason  et  al.,  2005;  Wilkinson  et  al.,  2005),  air  and  lichen 
samples  influenced  by  anthropogenic  Zn  (Cloquet  et  al.,  2006;  Dolgopolova  et  al.,  2006), 
cultured  phytoplankton  (Gelabert  et  al.,  2006),  and  seawater  (Bermin  et  al.,  2006). 

1.2  Zinc  stable  isotope  analysis 

Zinc  has  five  stable  isotopes:  64Zn  (49.2%),  66Zn  (27.8%),  h7Zn  (4.0%),  ftKZn  ( 1 8.4%),  and 
°Zn  (0.6%)  (Tanimizu  et  al.,  2002).  MC-ICP-MS  can  be  used  to  measure  the  ratio  of  Zn 
isotopes  in  a  sample.  We  measure  MZn,  6AZn,  and  f'sZn,  the  three  most  dominant  Zn 
isotopes.  All  isotopic  measurements  were  made  on  an  IsoProbe  MC-ICP-MS.  Typical 
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Figure  1.1.  Published  measurements  ofZn  isotopes  in  natural  samples  when  this  thesis 
began.  Data  is  from  (Marechal  et  aL,  2000),  except  for  data  on  microbial  sulfides  (  Archer 
and  Vance,  2002),  and  deep-sea  carbonates  (Pichat  et  aL,  2003). 


180*  210*  240'  270'  300'  330'  O'  30*  60*  00’  120“  150'  180 


Primary 

Productivity 


-60' 


ISO*  210*  240'  270*  300*  330'  0’  30'  60'  90*  120*  150’  ISO* 


0.5  0.9  1.3 


S66  Zn  %„ 

Figure  1.2.  The  Zn  isotope  composition  of  manganese  nodules  from  around  the  world. 
From  Marechal  et.  al.  (2000). 
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analytical  conditions  are  shown  in  Table  1. 


1.2.1  Conventions  for  reporting  Zn  isotope  ratios 

Zinc  isotope  data  is  typically  reported  in  terms  of  the  MZn/66Zn  ratio  of  the  sample 
in  relation  to  the  “Lyon  -  JMC  standard”  from  the  lab  of  Francis  Albarcde  (Marechal  et 
al.,  1999)  according  to  the  equation: 


'CZnrZn)^ 

\  Zn)jMCsuin<Uml 


\ 

1  *1000 

) 


Similarly,  the  ratio  of  6SZn/64Zn  is  reported  as: 


^ZnrZn)JMCM 


-1 


•  1000 


In  this  thesis,  isotope  effects  (or  any  offset  between  the  ftt,6Zn  of  two  samples)  are 
reported  as: 

AW’Z«  =  ?rZnA  -  8*2*1  B 

where  A  and  B  arc  the  two  samples.  This  terminology  is  potentially  confusing  because  A 
terminology  is  also  used  to  describe  non-mass  dependant  mass  fractionation  of  O,  S,  and 
Hg  isotopes,  e  notation,  often  used  to  describe  biological  isotope  effects  in  measurements 
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Typical  IsoProbe  operating  conditions: 

Sample  cone:  Aluminum 
Skimmer  cone:  Ni  “X-cone” 

Nebulizer:  MicroMist  glass  nebulizer 

Aspiration  rate:  75  pL/min,  free  draw 

Desolvator:  APEX  Q 

Desol vating  membrane:  none 


Cool  gas  flow:  13-14  L/min 

Intermediate  gas  flow:  0.9-1, 0  L/min 

Neb  gas  1 :  0.2-0.3  L/min 

Neb  gas  2:  0.7-0.9  L/min 

Sample  matrix:  2%  (vol.)  HN03 

Sample  uptake  time  before  measurement:  3  min 

Data  acquisition  time:  3  min 

Rinsing  time  between  measurements:  6  min 


Collector  Mass  Predominant  isotopes 


LI 

60 

60Ni 

Ax 

63 

6JCu 

HI 

64 

64Zn, 

H3 

65 

65Cu 

H4 

66 

6fiZn 

H5 

67 

67Zn 

H7 

68 

68Zn 

Table  1.1.  Typical  analytical  conditions  for  Zn  isotope  analysis  on  the  IsoProbe. 
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of  8i3C,  is  inappropriate  because  e  is  also  used  in  reporting  Nd  isotope  measurements  to 
signify  that  result  is  in  parts-per-ten  thousand.  ‘ A6  notation,  where: 

A  d^Zn  =  du'ZnA  -  d^ZnB 

was  first  used  for  C  isotopes,  and  is  attractive  because  it’s  meaning  is  intuitively  clear 
(Shacklcton  and  Pisias,  1984).  Unfortuantcly,  this  terminology  is  not  currently  used 
within  the  trace-metal  community  and  changing  common  notation  is  difficult.  Issues  of 
terminology  will  undoubdetly  become  more  complicated  if  mass  independent  isotope 
effects  arc  discovered  in  more  trace  metals  such  as  Fc. 


1.2.3  Mass-dependant  fractionation 

For  mass-dependent  fractionation  processes,  the  magnitude  of  6<  xZn  should  be 
approximately  twice  the  magnitude  of  6w’Zn  (see  Chapter  5  for  further  discussion).  The 
exact  relationship  between  these  two  quantities  depends  on  the  law  that  is  used  to 
describe  isotope  fractionation.  The  commonly  used  “exponential  law”  takes  the  form: 


\m2j 


where  the  ratio  of  two  isotopes  after  fractionation  (Rf)  is  related  to  the  initial  ratio  of  the 
two  isotopes  (R,)  as  a  function  of  the  isotopic  masses  (mi  and  m2)  and  a  fractionation 
factor  (P).  According  to  the  exponential  fractionation  law,  the  ratio  of  6f,sZn  to  b6,'Zn  will 
be  1.985. 
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1.2.4  Correcting  for  instrumental  mass  bias 

1.2.4. 1  Using  Cu  or  a  Zn  double-spike  to  correct  for  mass  bias 

During  the  analysis  of  Zn  isotopes  by  MC-ICP-MS,  large  isotopic  fractionations 
are  induced  in  the  spectrometer  between  sample  induction  and  isotope  signal 
measurement.  Several  different  schemes  arc  available  to  correct  for  this  mass  bias.  The 
simplest  involves  sample-standard  bracketing.  Assuming  that  the  magnitude  of 
instrumental  mass  bias  is  constant  between  samples  and  standards,  the  measured 
difference  in  the  isotope  ratios  will  be  equal  to  the  actual  difference  in  the  isotope  ratios. 

In  practice,  small  differences  in  the  matrix  between  samples  and  standards  often 
leads  to  changes  in  instrumental  mass  bias.  Instrumental  mass  bias  must  be  continuously 
monitored  in  order  to  correct  for  this  effect.  For  Zn,  both  Cu  normalization  and  a  double- 
spike  method  can  be  used.  Cu  has  only  two  stable  isotopes,  ^’Cu  (69.2%)  and  <,5Cu 
(30.8%).  By  spiking  both  samples  and  standards  with  Cu.  changes  in  instrumental  mass 
bias  can  be  monitored  by  looking  at  changes  in  the  measured  65CuA3Cu  ratio.  The 
exponential  mass  bias  fractionation  factor  (P)  calculated  for  Cu  cannot,  however,  be 
directly  applied  to  correct  for  Zn  fractionation  (Fig.  3).  Instead  an  empirical  relationship 
between  the  values  of  p  must  be  established  (Marechal  ct  ai„  1999). 

A  double-spike  can  also  be  used  to  correct  for  instrumental  mass  bias.  To  do  this, 
all  samples  arc  spiked  with  a  known  mixture  of  non-natural  Zn  isotopes,  typically  ^Zn 
and  Zn.  By  simultaneously  collecting  data  on  the  abundances  of  64Zn,  66Zn,  f,7Zn.  and 

6  S 

Zn  it  is  possible  to  calculate  both  the  concentration  ratio  of  spike  to  sample  and  the 
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mass  bias  of  the  instrument.  It  is  then  possible  to  calculate  the  isotope  composition  of  the 
original  sample. 

1. 2.4.2  The  comparative  advantages  of  Cu  and  Zn  double-spike  correction 

Correction  with  a  Cu  spike  and  a  Zn  double-spike  each  have  their  own  advantages 
and  disadvantages.  A  drawback  of  the  Cu  correction  scheme  is  that  Cu  and  Zn  behave 
differently  during  ICP-MS  analysis.  This  is  apparent  from  the  fact  that  values  of  fi  are 
different  for  these  two  elements.  While  the  Cu-correction  scheme  outlined  above  has 
been  rigorously  tested  and  demonstrated  to  work  well  for  standards  under  many 
conditions  (Archer  and  Vance,  2004;  Marechal  ct  al.}  1999),  it  cannot  be  proved  that  this 
correction  scheme  always  works  for  real  samples. 

Another  difference  is  that  the  Cu  correction  scheme  requires  complete  recovery  of 
Zn  from  the  sample.  Conversely,  a  Zn  double-spike  can  be  added  to  the  sample  before 
processing  to  account  for  fractionation  of  the  isotopes  during  sample  processing.  This 
simplifies  sample  processing,  especially  in  cases  such  as  the  analysis  of  Zn  isotopes  in 
seawater  where  achieving  complete  recovery  of  Zn  can  be  difficult. 

Finally,  the  correction  schemes  differ  in  that  only  Cu-correction  allows  a  check 
for  the  presence  of  isobaric  interferences  by  simultaneously  measuring  &'*Zn  and  bMZn. 
Isobaric  interferences  will  interfere  differently  with  the  measured  values  of  b^Zn  and 
&68Zn,  causing  them  to  fall  off  the  expected  linear  relationship  governed  by  mass- 
dependant  fractionation.  In  double-spike  analysis,  b^Zn  and  b(’*Zn  cannot  be  measured 
simultaneously,  so  isobaric  interferences  cannot  be  distinguished  from  changes  in  sample 
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o  Cu  JMC and 2n  JMC  - fzn*  feu 

•  Cu  NIST  976  and  Zn  JMC 


In  (f>*Zn/MZn) 

Figure  1.3.  During  the  course  of  a  day,  the  measured  isotope  ratios  change  as 
instrumental  mass  bias  (|i,  or  /  in  this  figure)changes.  However,  the  extent  of 
fractionation  of  these  two  elements  is  not  the  same  for  the  two  elements  as  predicted  by 
the  exponential  fractionation  law.  In  the  case  shown  here,  for  example,  the  exponential 
mass  bias  factor  for  Zn  varies  by  more  than  the  exponential  mass  bias  factor  for  Cu  (/Zn  > 
/cu).  An  empirically  determined  relationship  between  fzn  and  fcu  allows  one  to  make  a 
proper  correction  for  varying  instrumental  mass  bias,  (Marechal  et  al.,  1999). 


fractionation  or  instrumental  mass  bias.  Polyatomic  isobaric  interferences  may  occur 
when  sample  matrix  is  not  completely  eliminated.  Polyatomic  interferences  are  of  special 
concern  on  the  IsoProbe  where  an  Al  sampling  cone  (rather  than  Ni)  is  required  for 
proper  heat  dissipation.  An  2  Al4(lAr  interference  on  mass  67  is  present  with  typical 
signals  of  several  mV. 

During  data  collection  for  this  thesis,  a  Cu-correction  scheme  was  most  often 
used.  A  double-spike  was  used  to  measure  Zn  isotopes  in  six  hydrothermal  samples 
(Chapter  3)  and  several  standards.  Both  correction  schemes  produced  the  same  results, 
though  the  internal  error  for  the  Zn  double-spike  analysis  was  larger.  The  origin  of  this 
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error  is  not  known,  but  may  result  from  variations  in  the  magnitude  of27Al40Ar 
interferences. 

It  might  be  possible  to  achieve  the  advantages  of  the  double-spike  method, 
without  introducing  susceptibility  to  isobaric  interferences,  by  making  measurements  in 
medium  or  high-resolution  mode  where  interferences  can  be  separated  from  Zn  by  mass. 
This  requires  a  significantly  larger  samples  size  in  order  to  produce  the  same  signal 
intensity. 

1.3  Thesis  chapters 

Chapter  2  presents  data  on  the  Zn  isotope  ratios  found  in  many  anthropogenic 
materials.  Laboratory  standards  were  found  to  be  very  fractionated  (up  to  several  pemiil) 
from  typical  continental  material.  The  more  common  types  of  anthropogenic  Zn,  such  as 
pennies  and  health  products,  were  found  to  have  similar  values  of  6w'Zn  to  continental 
materials.  1  discuss  how  this  data  can  be  used  to  estimate  the  6(,<’Zn  of  many 
anthropogenic  Zn  sources  into  the  environment  and  oceans. 

Chapter  3  discusses  Zn  isotope  data  for  hydrothermal  vent  fluids  and  chimneys. 
Fluids  666Zn  values  were  measured  on  samples  from  several  different  hydrothermal  fields 
in  both  the  Pacific  and  Atlantic.  A  relationship  was  discovered  between  the  temperature 
of  hydrothermal  fluids  and  their  Zn  isotope  composition.  Cooler  fluids  (<250°C)  had 
566Zn  values  that  were  significantly  heavier  than  basalt  or  high-temperature  fluids  leading 
us  to  hypothesize  that  the  subsurface  precipitation  of  light  Zn  sulfides  was  associated 
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with  fluid  cooling.  Measurements  from  fluids  and  chimney  sulfides  taken  from  the  same 
vent  allow  us  to  measure  the  isotope  effect  for  hydrothermal  precipitation  in  situ. 

In  Chapter  4  we  present  the  results  from  diatom  cultures  grown  at  a  range  of  Zn 
concentrations  and  analyzed  for  Zn  isotopes.  By  carefully  distinguishing  between 
extracellular  Zn  and  intracellular  Zn,  we  have  been  able  to  measure  the  isotope  effects 
associated  with  separate  high-affinity  (A^Zn  =  -0.2%o)  and  low-affinity  (Ahf,Zn  -  -0.8%o) 
transport  pathways.  These  are  the  first  experiments  to  explain  transition  metal  isotope 
fractionation  in  terms  of  the  molecular  biology  of  uptake  mechanisms.  This  data  provides 
a  means  of  relating  the  &M,Zn  of  seawater  to  the  biological  cycling  of  Zn  in  the  oceans. 

Chapter  5  details  efforts  to  measure  Zn  isotopes  in  natural  seawater  and 
phytoplankton  samples.  I  have  developed  new  methods  for  the  analysis  of  Zn 
concentration  in  seawater,  and  a  new  method  for  the  Mg(OH)2  co-precipitation  of  Zn 
from  seawater  for  isotopic  analysis.  Analytical  problems  with  seawater  isotope  analysis 
arc  discussed.  6w’Zn  values  for  marine  plankton  tows  range  from  about  0.0%o  to  +0.6%o. 
Analysis  of  five  separate  samples  of  seawater  from  the  deep  Pacific  yielded  an  average 
&w'Zn  of+0.50±  0.08%o.  There  appears  to  be  a  trend  towards  lighter  Zn  values  of  &66Zn 
in  surface  seawater.  This  suggests  the  preferential  removal  of  isotopically  heavy  Zn  from 
the  surface  ocean.  I  hypothesize  that  adsorption  of  Zn  onto  sinking  particles  is  a  major 
pathway  for  the  removal  ofZn  from  the  surface  oceans. 

Chapter  6  briefly  explores  some  fundamental  questions  about  Zn  cycling  in  the 
ocean,  and  the  factors  that  control  Zn  isotope  distribution  in  the  oceans. 
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Chapter  2 


The  isotopic  composition  of  some  common 
forms  of  anthropogenic  zinc 


Anthropogenic  sources  account  for  much  of  the  zinc  (Zn)  in  the  environment.  Constraining  the 
isotopic  composition  of  anthropogenic  Zn  is  therefore  essential  to  understanding  the  environmental 
biogeochcmical  cycling  of  Zn  isotopes.  This  study  examines  the  isotopic  variability  in  several  different 
kinds  of  anthropogenic  Zn.  Laboratory  standards  are  highly  purified  and  can  be  significantly 
fractionated  from  natural  Zn.  Industrial  Zn  dust  and  U.S.  pennies,  which  arc  made  from  the  most 
common  grade  of  Zn  metal  (Special  High  Grade),  were  studied  to  represent  common  Zn  metal. 
Vitamins  were  studied  because  they  arc  more  highly  purified  than  typical  Zn  metals  and  arc  made  from 
chemical  compounds  such  as  Zn  oxide  or  Zn  gluconate.  The  isotopic  composition  of  laboratory 
standards  ranged  in  S66 Zn  from  -9. 1 5%o  to  +0. 1 7%o  compared  to  Lyon-JMC  Zn.  Zn  dust  and  pennies 
ranged  from  +0.09%o  to  +0.3  l%o,  and  vitamins  ranged  from  +0.09%o  to  +0.27%o.  External 
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rcproducibity  was  0.052%o  (2o  s.d.)  based  on  several  standards,  and  was  similar  to  internal  error.  The 
isotopic  range  for  non-laboratory  standards  is  much  smaller  than  the  total  range  seen  in  Zn  ore 
deposits,  but  has  a  similar  average  Zn  isotope  ratio.  This  data  presents  a  first  look  at  the  range  of  &MZn 
values  that  is  typical  of  common  anthropogenic  Zn  products. 

2.1  Introduction 

Zinc  (Zn)  is  commonly  used  in  many  man-made  items  and  over  3  million  metric  tons  arc 
estimated  to  be  released  into  the  environment  every  year  (Graedcl  ct  al.,  2005).  It  is  used  in  the 
manufacture  of  galvanized  steel,  in  alloys  with  other  metals  in  many  objects,  in  the  manufacture  of 
rubber  to  neutralize  acidity,  and  in  agriculture  as  a  crop  nutrient  (Gordon  et  al„  2003).  Zn  is  also  used 
in  many  health  products  and  is  a  common  component  in  vitamins,  sunscreen,  cold  medicine,  and  skin 
creams.  The  purification  ofZn  from  ores  as  described  by  (Gordon  et  al.,  2003)  is  a  multi-step  process. 
The  first  step  is  a  high  temperature  roast  to  convert  Zn  sulfides  into  oxides.  The  resulting  calcine  is 
then  dissolved  in  sulfuric  acid  and  impurities  may  be  removed  by  precipitation.  The  final  purification 
step  is  electrowinning,  in  which  high-purity  Zn  is  electroplated  from  solution  onto  aluminum  cathodes. 
Zn  may  then  be  manually  removed  from  the  cathodes  and  melted  for  processing  into  ingots,  shot,  dust, 
and  other  forms.  A  less  common  method  of  Zn  refining  uses  high-temperature  distillation  of  metallic 
Zn  in  place  of  electrowinning.  Zn  oxide  is  typically  made  by  evaporation  and  oxidation  of  Zn  and  Zn 
oxide  is  a  precursor  to  many  other  Zn  compounds  (FDA,  2002), 

Anthropogenic  sources  arc  a  major  contributor  to  the  Zn  found  in  natural  environments. 
Approximately  half  of  total  Zn  emissions  to  the  atmosphere  are  thought  to  have  anthropogenic  sources 
(Pacyna  and  Pacyna,  2001 )  and  as  much  as  half  of  the  global  fluvial  flux  of  Zn  to  the  oceans  may  be 
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anthropogenic  (Shiller  and  Boyle,  1985).  The  influence  of  anthropogenic  Zn  is  seen  even  in  remote 
locations  far  from  the  major  sources  of  Zn  pollution.  For  example,  one  third  of  the  Zn  dry  deposition  to 
the  ocean  near  Bermuda  is  thought  to  be  anthropogenic  (Arimoto  et  al.,  2003),  and  anthropogenic  Zn 
has  been  measured  in  ice  from  both  Greenland  (Candelone  et  al.,  1995)  and  Antarctica  (Planchon  ct  al., 
2002).  In  order  to  model  the  global  cycling  of  Zn  isotopes,  wc  need  to  characterize  the  isotopic 
composition  of  anthropogenic  Zn. 

Previous  studies  have  found  large  variations  in  the  Zn  isotope  composition  of  anthropogenic  Zn 
standards  (Mason  et  al.,  2004;  Tanimizu  et  al.,  2002).  The  standards  measured  by  Tanimizu  et  al. 
ranged  from  6M'ZnjMc  =  -2.4 1  %o  to  +O.I3%o.  This  suggests  that  isotopic  fractionation  may  be  common 
during  the  purification  of  laboratory  standards.  Tanimizu  ct  al.  found  that  the  sample  with  the  greatest 
deviation  from  continental  material  was  NIST-SRM  682,  which  was  also  the  purest  Zn  standard  tested. 
They  suggest  that  this  large  fractionation  resulted  from  the  distillation,  zone  refining,  and  gasification 
employed  in  the  final  purification  of  this  sample.  Metal  reduction  during  electroplating  may  also  have 
a  significant  kinetic  isotope  effect  (Kavner  et  al.,  2005),  suggesting  that  electrolytically-purified  Zn 
could  also  have  a  unique  anthropogenic  signature. 

N on-electrolytic  processes  may  also  fractionate  Zn  isotopes.  Analysis  of  Zn  isotopes  samples 
from  a  Pb-Zn  refinery  and  a  steel  mill  show  variations  in  6w’Zn  from  -0.63%o  to  +0.58%o.  In  the  Pb-Zn 
refinery,  there  was  a  trend  towards  lighter  isotope  values  near  the  end  of  the  refining  process,  with  the 
lightest  value  found  in  the  chimneystack  (Matticlli  et  al.,  2005).  Zn  isotopes  from  lichens  collected 
near  an  ore-processing  and  mining  site  arc  heavier  than  minerals  collected  from  the  mining  site  or  the 
typical  mineralogical  spread  of  Zn  isotope  values  (Dolgopolova  ct  al.,  2006).  The  authors  suggest  that 
such  isotopic  signals  may  be  used  to  trace  anthropogenic  or  natural  Zn  in  aerosols.  Zn  isotopes  in 
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rainwater  in  France  varied  by  about  0.3%o,  and  differences  were  seen  between  anthropogcnically 
influenced  and  marine-source  rain  (Luck  et  al.,  1999). 

Mineral  ores  may  also  be  a  source  of  variability  in  anthropogenic  Zn.  Samples  collected  at 
ancient  hydrothermal  sites,  the  source  of  modern  Zn  ores,  range  from  6fif’Zn  -  -0.43%o  to  +  1.33%o 
(Mason  et  al.,  2005;  Wilkinson  et  al.,  2005).  Wilkinson  et  al.  compiled  Zn  isotope  data  for  all 
terrestrial  samples  measured  to  date  to  show  that  they  span  the  entire  range  between  the  heaviest  and 
lightest  hydrothermal  sample.  Three  homogenized  ore  samples  have  a  smaller  isotopic  range,  between 
6w’Zn  =  -0.06%o  and  +0.33%o  (Chapman  et  al.,  2006).  Zn  isotopes  arc  variable  on  spatial  scales  of 
centimeters  and  millimeters  (Mason  et  al.,  2005),  so  the  lesser  variability  found  in  ores  may  represent 
the  effects  of  homogenization. 

While  the  total  variability  in  anthropogenic  Zn  isotopes  including  laboratory  standards  and 
aerosols  may  be  several  permil,  the  variability  ofZn  isotopes  in  the  most  common  forms  of 
anthropogenic  Zn  has  not  previously  been  investigated.  Several  grades  ofZn  arc  manufactured  and 
used  in  the  US  including  Special  High  Grade  (SHG),  High  Grade  (HG),  Prime  Western  (PW)  with 
minimum  purities  of  99.99%,  99.95%  and  98.50%,  respectively.  Together,  these  three  grades  account 
for  94%  of  the  reported  Zn  use  in  the  US  during  2004  and  SHG  Zn  alone  accounted  for  63%  of  Zn  use 
during  this  time  (Jorgenson,  2004).  U.S.  pennies,  which  have  been  made  from  Zn  with  a  thin  coating  of 
Cu  since  1982,  arc  made  with  SHG  Zn  and  arc  therefore  representative  of  the  majority  of  Zn  used  in 
the  US  (Jasinski,  1994).  If  there  is  isotopic  fractionation  associated  with  industrial  Zn  refining  and 
production  processes,  we  expect  this  to  be  most  apparent  in  SHG  Zn  compared  to  less  purified  HG  and 
PW  grade  Zn.  We  have  also  measured  some  health  and  medical  products;  although  they  represent  only 
a  small  portion  ofZn  use.  they  require  the  most  highly  purified  Zn  used  in  common  anthropogenic 


32 


products  and  are  made  from  Zn  compounds  such  as  Zn  oxide  and  Zn  gluconate  rather  than  Zn  metal 
(FDA,  2002). 

2.2  Materials  and  methods 

2.2.1  Zinc  samples 

Dissolved  Zn  standards  were  N1ST-SRM  682  and  a  1000  ppm  Zn  standard  (AccuTracc, 
AccuStandard  Inc.,  Lot#B2075078).  Laboratory-grade  Zn  metal  samples  were  10  mesh  (#1:  B&A,  Lot 
D355Z077R)  and  30  mesh  (#2:  J.T.  Baker,  Lot  G42701).  Laboratory  Zn  acetate  standards  were  from 
Mallinckrodt  (#1:  Lot  WALP,  and  #2:  Lot  KDSJ)  and  Baker  Scientific  (#3:  Lot  H20157).  Industrial  Zn 
metal  dust  samples  of  98.5%  purity  were  obtained  from  the  US  Zinc  (Austin,  TX)  as  regular  (USZ  # 

I),  regular  superfine  (USZ  #  5),  and  low  Pb  Zn  (USZ  #1  XL).  Industrial  high  purity  Zn  shot  (99.995%) 
was  obtained  from  the  Canadian  Electrolytic  Zinc  (CEZ)  refinery  (Salaberry-dc-Vallcyficld,  Quebec) 
through  Falconbridge  Limited.  U.S.  pennies  were  obtained  from  circulation  and  arc  noted  by  the  year 
in  which  they  were  minted.  Health  and  medical  products  were  obtained  from  local  commercial  sources 
and  include  Benadryl  Itch  Stopping  Cream  (Pfizer  Consumer  Healthcare,  Lot  #23454L),  Centrum 
Silver  multivitamins  (Wyeth  Consumer  Healthcare,  Control  #A98404),  Dr.  Zinc  lozenges  (McKesson, 
Lot  #3KVO305),  Nature’s  Bounty  zinc  tablets  (Nature’s  Bounty,  Inc.,  Lot  #20515101),  and  ColdEczc 
zinc  lozenges  (The  Quigley  Corporation). 

2.2.2  Sample  preparation 

Samples  were  handled  under  Class- 100  clean  air  flow  conditions.  Water  was  deionized  and 
distilled  in  a  borosilicatc  glass  still  (Corning  McgaPure),  and  all  acids  were  triply-distilled  in  Vycor 
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glass  to  remove  contaminants.  All  labwarc  was  acid  cleaned  PFA  (Savillex),  polypropylene  or 
polyethylene.  Metal  samples  were  prepared  by  dissolving  directly  in  2%  HNO.i  for  analysis,  or  by 
dissolving  in  2N  HCl  prior  to  purification  by  anion  exchange  chromatography.  Zn  acetate  samples 
were  twice  dissolved  in  2%  HNO3  and  evaporated  to  drive  off  acetic  acid  before  redissolution  in  2% 
HNOi.  Pennies  were  washed  with  2%  Citranox  detergent,  rinsed  with  distilled  water,  and  dissolved  in 
2N  HCl  prior  to  anion  exchange  purification.  Health  and  medical  products  were  dissolved  in  2%  nitric 
acid  and  small  aliquots  were  removed  and  evaporated  to  dryness.  These  subsamples  were  then 
combusted  for  8  hours  at  450,,C  to  remove  all  organic  material  and  the  remaining  material  was 
dissolved  in  2N  HCl  prior  to  anion  exchange  purification. 

Between  10  to  200  pg  of  each  sample  was  purified  by  anion  exchange  chromatography  on  AG- 
MP1  resin  (BioRad)  according  to  a  method  previously  used  for  the  purification  of  Fe,  Cu,  and  Zn  for 
isotopic  analysis  (Marcchal  ct  al.,  1999).  Because  samples  were  not  being  processed  for  Fc  and  Cu 
isotope  analysis,  samples  were  loaded  in  2N  HCl.  The  sample  was  rinsed  with  15  mL  of  2N  HCl  to 
elute  other  elements,  and  then  the  Zn  fraction  was  eluted  in  1 2  mL  0.5N  I  INOi.  The  eluent  was 
evaporated  and  reacted  with  200  pL  of  14N  HNOj  and  100  pL  of  H2O2  (Ultrcx  II,  J.T.  Baker),  to 
oxidize  organics  that  may  have  leached  from  the  resin.  With  every  group  of  samples  processed  by 
column  chromatography,  two  process  blanks  and  one  or  two  isotope  standards  were  processed  and 
analyzed  to  assess  the  magnitude  of  contamination  and  matrix  effects  during  measurement. 

2.2.3  Isotope  analysis 

All  samples  were  diluted  to  200  ppb  Zn  in  2%  HNO3  for  isotope  analysis  and  spiked  with  100 
ppb  Cu  (Ultra  Scientific,  Lot  #D00204)  to  monitor  instrumental  mass  bias.  Data  was  collected  on  an 
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IsoProbc  multi-collector  ICP-MS  (Thermo  Electron,  formerly  Micromass),  using  Faraday  collectors  to 
monitor  signals  on  masses  60,  63,  64,  65,  66,  67,  and  68.  Samples  were  introduced  into  the  plasma  by 
an  Apex  Q  inlet  system  with  a  dcsolvating  membrane  (ESI)  using  a  75  pL/min  MicroMist®  (Glass 
Expansion)  nebulizer.  On-peak  acid  blank  subtraction  was  applied  to  all  samples  by  monitoring  the 
signal  from  2%  HNO?  during  one  60  second  block,  three  minutes  before  sample  analysis.  Sample  and 
standard  data  was  collected  in  fifteen  cycles  of  ten  seconds,  and  up  to  two  cycles  more  than  two 
standard  deviations  from  the  mean  were  discarded. 

The  only  spectral  interference  for  which  we  apply  a  correction  is  64Ni  (overlapping  with  hJZn). 
This  correction  is  made  by  monitoring  the  abundance  of  w,Ni,  assuming  natural  isotopic  abundances, 
and  correcting  for  instrumental  mass  bias  with  an  internal  Cu  spike  (sec  below).  Ni  corrections  were 
insignificant  at  our  analytical  precision.  We  first  determine  Zn  isotope  ratios  by  sample-standard 
bracketing  alone  and  then  correct  these  data  for  sample  matrix  effects  on  instrumental  mass  bias  using 
an  internal  Cu  spike.  The  relationship  between  Zn  and  Cu  mass  bias  can  differ  from  that  expected 
based  on  the  exponential  mass  bias  law  alone  (Albarcdc  and  Beard,  2004;  Marcchal  et  al.,  1999). 
Following  their  suggestion,  we  determine  the  linear  relationship  between  the  natural  log  of  measured 
^ZnA'Zn  versus  measured  ft5CuA3Cu  for  all  standards  analyzed  during  a  single  session,  a  procedure 
termed  empirical  external  normalization  (EEN)  (Clayton  et  al.,  2002;  Mason  ct  al.,  2004).  We  apply 
this  correction  to  both  samples  and  standards,  for  a  'modified'  EEN  sample-standard  bracketing 
technique.  All  samples  were  measured  in  triplicate  and  ^Zn/^Zn  ratios  were  measured  to  verify  mass- 
dependent  fractionation  and  the  absence  of  Zn  spectral  interferences. 
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NIST-SRM  682  has  been  measured  repeatedly  against  the  Lyon-JMC  Zn  (Marcchal  ct  a!., 

1999)  and  is  lighter  by  866ZnjMc  =  -2.45  %o.  All  measurements  arc  corrected  with  this  offset  so  that 
isotope  ratios  may  be  reported  compared  to  Lyon-JMC  Zn. 

2.3  Results  and  Discussion 
2.3.1  Error  analysis 

The  measured  isotope  ratio  of  samples  compared  to  unprocessed  bracketing  standards  can  be  affected 
by  Zn  contamination  of  the  samples,  isotopic  fractionation  during  the  purification  procedure,  or 
analytical  changes  resulting  from  matrix  residues.  We  discount  the  effect  of  Zn  contamination  because 
process  blank  concentrations  were  low  (consistently  between  0  and  10  ng  and  never  more  than  0. 1%  of 
the  sample  concentrations).  Assuming  that  contamination  sources  have  an  isotopic  ratio  similar  to 
other  continental  and  anthropogenic  samples,  this  amount  of  contaminant  Zn  is  too  small  to  have  a 
measurable  effect  on  the  isotope  ratio  of  samples.  Zn  recoveries  were  always  greater  than  98%,  so 
there  should  be  no  isotopic  fractionation  associated  with  column  purification  and  other  processing 
steps.  NIST-SRM  682  process  standards  were  treated  by  the  same  purification  protocol  as  samples  in 
order  to  gauge  the  accuracy  and  precision  of  our  analytical  methods  (Figure  1 ).  Most  of  the  Cu- 
corrected  ratios  are  higher  than  the  raw  data  corrected  by  sample-standard  bracketing  alone.  Although 
the  magnitude  of  this  correction  is  small  (the  average  A86f,Zn  =  +0,02%o  for  the  five  standards),  it  may 
reflect  a  genuine  systematic  matrix  effect  on  the  Cu  isotope  ratio  in  samples. 

In  our  experience,  though,  the  main  source  of  analytical  error  is  occasional  appearance  of  a 
systematic  offset  between  pre-column  and  post-column  standards,  as  has  also  been  observed  in  other 
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Figure  2.1.  Sample-standard  bracketing  corrected  (O)  and  sample  corrected  (•)  Zn  isotope  ratios  of 
standards  processed  alongside  samples  on  several  different  dates. 


studies  (Archer  and  Vance,  2004).  In  some  analytical  sessions,  use  of  the  dcsolvating  membrane 
attached  to  the  APEX  resulted  in  systematic  offsets  in  all  post-column  samples  of  as  much  as  +0.2%o. 

In  this  study,  careful  post-column  treatment  of  the  samples  with  HNO3/H2O2  to  remove  organics  seems 
to  have  largely  prevented  this  problem.  In  other  experiments,  we  found  that  removing  the  dcsolvating 
membrane  appears  to  eliminate  this  problem. 

Reproducibility  for  triplicate  analysis  of  individual  samples  was  between  0.00%o  and  0.06%o 
(2 o  s.d.)  (Table  1).  The  difference  in  measured  isotope  ratio  (A)  between  standards  used  for  sample- 
standard  bracketing  and  standards  run  through  the  full  sample  purification  protocol  was  used  to 
estimate  external  reproducibility.  The  standard  deviation  of  Aw’Zn  for  the  five  process  standards  was 
0.052%o  (2o),  greater  than  the  internal  reproducibility  within  triplicate  analyses  for  nearly  all  of  our 
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Figure  2.2.  Zn  isotopic  composition  of  anthropogenic  samples  showing  A)  the  entire  measured  range 
and  B)  the  smaller  range  in  which  all  non- laboratory  standards  fall. 


samples.  We  are  wary  of  applying  the  error  calculated  for  post-column  standards  to  samples  with  a 
potentially  more  complicated  matrix,  but  we  believe  that  this  error  is  more  appropriate  than  the  smaller 
errors  calculated  from  internal  reproducibility.  A  standard  deviation  of0.052%o  is  equivalent  to  a 
standard  error  of  0.030%o  (2a)  for  triplicate  analysis.  Several  samples  were  run  in  different  analytical 
sessions  months  apart,  these  replicate  analyses  give  the  same  S^Zn  values  within  the  analytical  error. 

2.3.2  Zn  isotopes  in  anthropogenic  samples 

Zn  isotope  ratios  were  measured  for  twenty-two  anthropogenic  Zn  samples  (Figure  2).  The  isotopic 
compositions  of  the  laboratory  standards  ranged  from  b^ZnjMr  =  -9. 15%o  to  +0. 1 7%o.  The  AccuTracc 
standard  has  the  lightest  Zn  isotope  ratio  reported  to  date.  N1ST-SRM  682  and  two  of  the 
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Sample 

name 

Sample-Standard 
Bracketing  (SSB) 

corrected  data 

Sample-Standard  Bracketing/ 

Cu  Empirical  External  Normalization 
(EEN)  corrected  data 

S^Zn^sT  666 

-SRM  682  JMU 

s.d.  (2a) 

&66ZnN|ST 

-SRM  662 

&66ZnJMC 

s.d.  (2a) 

NIST  SRM-682 

0.00 

-2.45 

AccuTrace  Std 

-6.70 

-9.15 

Metal  (10  mesh 

2.61 

0.16 

0.03 

2.62 

0.17 

0.00 

Metal  (30  mesh 

2.53 

0.08 

0.03 

2.54 

0.09 

0.03 

Acetate  #1 

2.61 

0.16 

0.24 

2.46 

0.01 

0.04 

Acetate  #2 

2.61 

0.16 

0.10 

2.53 

0.08 

0.05 

Acetate  #3 

2.51 

0.06 

0.06 

2.39 

-0.06 

0.04 

Std  Dust 

2.41 

-0.04 

0.09 

2.60 

0.15 

0.02 

Fine  Dust 

2.46 

0.01 

0.02 

2.64 

0.19 

0.04 

Low  Pb  dust 

2.48 

0.03 

0.08 

2.54 

0.09 

0.03 

CEZ  shot 

2.66 

0.21 

0.02 

2.67 

0.22 

0.03 

1 983  Penny 

2.51 

0.06 

0.07 

2.60 

0.15 

0.02 

1 983  Penny  #2 

2.56 

0.11 

0.03 

2.59 

0.14 

0.04 

1985  Penny 

2.54 

0.09 

0.07 

2.65 

0.20 

0.04 

1985  Penny  #2 

2.56 

0.11 

0.03 

2.63 

0.18 

0.03 

1991  penny 

2.49 

0.04 

0.11 

2.67 

0.22 

0.06 

1991  Penny  #2 

2.59 

0.14 

0.03 

2.64 

0.19 

0.05 

1999  Penny 

2.55 

0.10 

0.06 

2.62 

0.17 

0.01 

2000  Penny 

2.69 

0.24 

0.07 

2.76 

0.31 

0.05 

2000  Penny  #2 

2.69 

0.24 

0.02 

2.73 

0.28 

0.02 

2002  Penny 

2.61 

0.16 

0.06 

2.65 

0.20 

0.05 

Benadryl 

2.54 

0.09 

0.01 

2.56 

0.11 

0.02 

Centrum 

2.53 

0.08 

0.05 

2.54 

0.09 

0.02 

Dr  Zn 

2.63 

0.18 

0.03 

2.65 

0.20 

0.02 

Nature's  Bounty 

2.59 

0.14 

0.24 

2.55 

0.10 

0.03 

ColdEeze 

2.61 

0.16 

0.07 

2.72 

0.27 

0.04 

Tabic  2.1.  Zn  isotopic  composition  of  anthropogenic  samples  compared  to  NIST  SRM-682  and  Lyon- 
JMC  Zn.  Standard  deviations  arc  reported  for  triplicate  analysis  of  individual  samples. 
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Zn  acetate  samples  (#1  and  #3)  were  lighter  than  any  of  the  non-laboratory  Zn  samples  measured.  The 
S^’Zn  of  Zn  dust  ranged  from  +0.09%o  to  +0.22%o,  6f  <’Zn  for  pennies  ranged  from  0.1 4%o  to  0.3  l%o, 
and  the  health  and  medical  products  566Zn  ranged  from  0.09%o  to  0.27%o. 

All  common  types  of  anthropogenic  Zn  (excluding  laboratory  standards)  fell  within  a  small 
range  of  isotope  compositions.  The  total  range  in  8ftf,Zn  that  we  measured  for  common  Zn  is  +0.09%<> 
to  +0.28%o,  more  than  an  order  of  magnitude  smaller  than  the  range  seen  in  hydrothermal  samples 
discussed  in  Section  1.  Interestingly,  the  average  isotopic  composition  of  our  “common”  Zn  is  8(rf,Zn  - 
+0.1 8%o,  quite  similar  to  the  +0.l5%o  average  isotope  composition  of  samples  measured  in  two  Zn  ore 
fields  (Mason  et  al.,  2005;  Wilkinson  ct  al.,  2005).  Based  on  this,  we  suggest  that  the  S6hZn  variability 
we  see  in  different  anthropogenic  samples  represents  isotopic  differences  between  ores,  modified  by 
homogenization  during  processing.  We  do  not  find  evidence  of  a  large  isotope  effect  associated  with 
the  purification  ofZn  ore  into  metal. 

2.4  Conclusions 

When  studying  Zn  isotopes  in  the  environment  (for  example,  in  seawater  (Bcrmin  et  al.,  2006), 
marine  materials  (Marechal  et  al.,  2000),  or  plants  (Weiss  et  al.,  2005),  there  may  be  both  natural  and 
anthropogenic  sources  of  Zn.  We  believe  that  the  samples  analyzed  here  are  broadly  representative  of 
both  the  average  isotopic  composition  and  the  variability  that  will  be  found  in  common  forms  of 
anthropogenic  Zn. 
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Chapter  3 


Zn  isotopes  in  hydrothermal  vent 
fluids  and  chimneys 


Many  of  the  heaviest  and  lightest  natural  zinc  (Zn)  isotope  ratios  have  been 
discovered  in  hydrothermal  ore  deposits.  However,  the  processes  responsible  for 
fractionating  Zn  isotopes  in  hydrothermal  systems  arc  poorly  understood.  We  have 
studied  Zn  isotopes  in  seafloor  hydrothermal  chimneys  and  vent  fluids  from  several  vent 
fields  along  the  East  Pacific  Rise  (EPR)  and  Mid-Atlantic  Ridge  in  order  to  better  assess 
the  total  range  of  Zn  isotopes  in  hydrothermal  systems  and  to  understand  the  factors 
which  arc  responsible  for  this  isotopic  fractionation.  Zn  isotopes  may  then  be  used  to 
better  understand  hydrothermal  processes.  Of  all  the  physical  and  chemical  parameters 
examined  such  as  pH  and  concentrations  of  many  trace  elements,  only  temperature  was 
found  to  correlate  with  fluid  6w,Zn  values.  Lower  temperature  fluids  (<250C)  had  both 
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heavier  and  more  variable  6M,Zn  values  compared  to  higher  temperature  fluids  from  the 
same  hydrothermal  field.  Similarly,  chimney  sulfides  from  a  low-temperature  vent  at 
EPR  9°N  were  found  to  be  heavier  and  more  variable  in  66f>Zn  values  than  sulfides 
sampled  from  a  high-temperature  vent.  We  hypothesize  that  subsurface  cooling  of 
hydrothermal  fluids  leads  to  the  precipitation  of  isotopically  light  Zn  in  sulfides,  and  that 
this  process  is  the  primary  cause  of  Zn  isotope  variation  in  hydrothermal  fluids. 

3.1  Introduction 

With  multi-collector  inductively  coupled  plasma  mass  spectrometry  (ICPMS), 
stable  isotope  variability  of  transition  metals  has  been  seen  in  many  natural  samples. 
Metal  isotopes  studied  in  active  hydrothermal  systems  include  Fe  isotopes  (Rouxcl  et  al., 
2004b;  Severmann  ct  al.,  2004;  Sharma  ct  al.,  2001),  Cu  isotopes  (Rouxcl  ct  al.,  2004a; 
Zhu  et  al.,  2000),  Mo  isotopes  (McManus  et  al.,  2002),  Se  isotopes  (Rouxcl  ct  al.,  2004b), 
and  Sb  isotopes  (Rouxcl  ct  al.,  2003). 

Zn  isotopes  have  been  studied  in  rocks  (Chapman  et  al.,  2006;  Marechal  et  al., 
2000),  marine  sediments  (Pichat  et  al.,  2003),  biological  materials  (Marechal  ct  al.,  2000; 
Weiss  et  al.,  2005),  seawater  (Bcrmin  ct  al.,  2006),  and  ore  deposits  (Mason  ct  al.,  2005; 
Wilkinson  et  al.,  2005).  It  is  possible  that  Zn  isotopes  may  help  contribute  to 
understanding  ore  formations  and  hydrothermal  processes.  Economic  deposits  ofZn  ore 
originate  mainly  as  hydrothermal  vent  fields  (Herzig,  1999),  and  Zn  isotopes  have  been 
proposed  as  a  tool  for  distinguishing  between  Zn-rich  areas  in  mining  deposits  and  distal 
regions  where  Zn  is  less  concentrated  (Wilkinson  et  al.,  2005). 
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The  fractionation  of  Zn  isotopes  in  hydrothermal  systems  may  influence  the  Zn 
isotope  ratios  seen  in  the  environment.  The  deposits  from  hydrothermal  systems  contain 
both  the  heaviest  and  lightest  Zn  isotope  samples  reported  to  date  (6<>6Zn  =  -0.43%o  and 
6h6Zn  =  +1.33%o)  (Wilkinson  et  al.,  2005).  Zn  isotope  signals  may  be  transported 
throughout  the  environment  by  the  weathering  of  ancient  hydrothermal  deposits  and 
contemporary  hydrothermal  fluid  discharge. 

Thorough  studies  of  Zn  isotopes  have  been  carried  out  in  the  Alexandrika 
volcanic-hostcd  massive  sulfide  ore  deposit  in  Urals,  Russia  (Wilkinson  et  al.,  2005)  and 
the  Irish  Midlands  ore  field  (Mason  et  at.,  2005).  The  Alexandrika  ore  deposit  formed  as 
a  seafloor  hydrothermal  system  similar  to  modern  systems  studied  here.  Samples  were 
analyzed  from  the  hydrothermai-metasomatic  stockwork  believed  to  be  the  feeder  zone  to 
the  hydrothermal  deposit,  a  seafloor  chimney,  and  some  clastic  sediments.  66ftZn  values 
from  the  deposit  ranged  from  -0.43%o  to  +0.23%o.  A  systematic  increase  in  6hf,Zn  values 
from  the  chimney  core  to  the  chimney  rim  was  attributed  to  either  a  temperature 
dependence  in  the  fractionation  factor  for  precipitation  or  Rayleigh  distillation  as  fluids 
diffuse  through  the  chimney  wall.  Zn  isotopes  in  the  Irish  Midlands  ore  field  show  a 
general  trend  of  lighter  Zn  isotopes  precipitated  in  the  deep  feeder  veins  and  heavier  Zn 
isotopes  near  the  top  of  the  hydrothermal  system.  This  data  is  consistent  with  Rayleigh 
distillation  of  Zn  with  isotopically  light  Zn  partitioning  into  Zn  sulfides.  Other  studies 
have  investigated  Zn  isotopes  in  sulfides  from  various  ore  deposits  (Marechal  et  al., 

1999)  and  Archcan  shales  (Archer  and  Vance,  2002). 
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Despite  these  initial  studies,  the  cause  ofZn  isotope  fractionation  in  hydrothermal 
systems  remains  poorly  constrained.  It  is  hoped  that  the  study  of  active  hydrothermal 
vents  will  permit  to  provide  a  better  understanding  of  how  factors  such  as  fluid 
temperature,  temperature  of  sulfide  precipitation,  and  the  mineralogy  of  hydrothermal 
precipitates  will  affect  Zn  isotope  fractionation. 

3.1.1  Geologic  settings 

The  21°N  site  on  the  East  Pacific  Rise  (EPR)  was  among  the  first  ridge-crest 
hydrothermal  systems  to  be  studied  (Von  Damm  ct  al,  1 985).  Three  of  the  vent  sites 
(Ocean  Bottom  Seismometer  (OBS),  South  West  (SW),  and  Hanging  Gardens  (HG)) 
were  typical  high  temperature  black  smokers  and  their  chemical  composition  remained 
quite  stable  over  the  three  years  between  1979  and  1 981  (Campbell  et  al.,  1988).  The 
fourth  (National  Geographic  Society  (NGS))  was  a  black  smoker  when  originally 
sampled  in  1979.  The  vent  had  become  scaled  off  when  revisited  in  1981,  and  it  emitted 
low-temperature  white  smoker  fluids  when  broken  open.  The  9DN  site  is  one  of  the 
youngest  hydrothermal  sites  on  the  EPR  as  determined  by  volcanic  activity  and 
-|<Jp0/-M,pb  dating  (Von  Damm.  2000).  The  site  was  first  sampled  in  1991  and  with  only 
fifteen  years  between  the  first  hydrothermal  activity  and  sampling,  wc  expected  less 
hydrothermal  alteration  here  than  at  other  sites.  The  Guaymas  basin  lies  in  the  Gulf  of 
California,  also  on  the  East  Pacific  Rise  but.  in  contrast  to  the  9°N  and  21°N  sites,  it  is 
covered  by  a  thick  layer  of  sediments  through  which  hydrothermal  fluids  pass  before 
venting  on  the  seafloor.  The  TAG  hydrothermal  field  lies  on  the  more  slowly  spreading 
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Mid- Atlantic  Ridge.  This  site  is  distinguished  by  the  large  extent  of  subsurface  Zn  sulfide 
redissolution,  particularly  in  the  numerous  white  smokers. 

3.2  Methods 

3.2.1  Sample  collection 

Most  hydrothermal  vent  fluids  were  obtained  from  the  collection  of  John  Edmond 
stored  at  the  Massachusetts  Institute  of  Technology  (MIT).  All  these  samples  were 
collected  using  a  titanium  syringe  with  a  snorkel  that  could  reach  directly  into  the  vent 
chimney  for  sampling  (Von  Damm  et  al.,  1985).  It  was  believed  at  the  time  that  this 
sampling  method  prevented  the  precipitation  of  sulfides  from  the  samples,  and  that  when 
sulfides  did  form  they  could  be  quickly  redissolved  by  the  addition  of  6N  HCI  to  a  pH  of 
about  1 .6.  In  fact,  these  early  sampling  methods  may  have  left  some  particles  behind  in 
the  syringes  (“dregs").  All  samples  were  visually  inspected  to  ensure  that  no  visible 
precipitates  remained.  Precipitates  that  formed  since  sample  collection  were  analyzed  to 
ensure  that  they  contained  very  little  Zn. 

For  fluids  from  the  MIT  hydrothermal  vent  collection,  the  highest  recorded 
temperature  from  the  vent  field  was  used  instead  of  the  temperature  recorded  for  any 
individual  sample.  Individual  temperature  measurements  were  quite  difficult  with  the 
ALVIN  temperature  probe  (Campbell  et  al.,  1988).  For  samples  collected  at  the  21°N  site, 
only  sample  1 160-7  had  a  recorded  temperature  (224°C)  that  was  more  than  a  few 
degrees  different  from  the  field  end-member  temperature  (350°C).  This  particular  sample 
most  likely  represents  an  inaccurate  temperature  measurement  (K.  Von  Damm,  personal 
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communication).  Trace-mctal  concentrations  in  this  sample  (including  Cu)  support  our 
assumption  that  this  is  a  high-temperature  vent. 

Samples  from  the  Pacific  were  collected  from  ALVIN  dives  during  1981-2  cruise 
to  the  hydrothermal  fields  at  2 1  °N  on  the  East  Pacific  Rise  and  the  Guay  mas  basin. 
Numerous  geochemical  parameters  were  reported  in  previous  studies,  including 
temperature,  pl-1,  alkalinity,  and  the  concentrations  of  Mg,  SiC>2,  Li,  Na,  K,  Rb,  Sc,  Ca, 

Sr,  Ba,  Cl,  SiC>4,  HjS.  NH*.  NOi,  PO4,  and  the  trace  metals  Mn,  Fc,  Co,  Cu,  Zn,  Ag,  Cd, 
Pb.  As,  Se,  and  A1  (Von  Damm,  1984). 

Samples  from  the  Mid- Atlantic  Ridge  were  collected  at  the  TAG  and  Snake  Pit 
hydrothermal  sites  in  1990.  Vents  sampled  on  this  cruise  included  both  high  temperature 
“black  smokers”  and  lower  temperature  “white  smokers".  White  smoker  and  black 
smoker  end-member  concentrations  (Edmond  ct  al.,  1991)  and  data  on  dissolved  trace- 
element  fluxes  have  been  previously  reported  (German  ct  al..  1991).  Chemical  and 
physical  data  on  individual  samples  was  provided  by  Andy  Campbell  (personal 
communication). 

Samples  from  the  hydrothermal  field  at  9-1 0°N  on  the  East  Pacific  Rise  were 
collected  in  2004  in  the  ALVIN  submersible  dives  4053  and  4057  during  the  cruise 
ATI  1-2  (Rouxcl  ct  al.,  in  prep).  The  high-temperature  vent  Bio9”  is  located  at  9‘’50'N 
within  an  area  where  dramatic  changes  in  the  hydrothermal  systems  have  taken  place. 
Prior  to  1994,  Bio9  was  the  only  high  temperature  vent  sampled  in  the  Bio9  area.  By 
1994,  the  Bio9’  vent  was  a  distinct  individual  black  smoker.  By  1999,  the  number  of 
smokers  around  Bio9  had  begun  to  increase,  replacing  some  of  the  areas  that  had 
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previously  been  characterized  by  diffuse  flow.  Bio9"  is  the  third  distinct  smoker  in  the 
Bio9  area  which  has  been  initially  sampled  in  2002.  K-vcnt  at  9’29.8°N  is  the  most 
southerly  active  vent,  before  encountering  a  zone  of  extinct  vents  (Von  Damm,  2000).  K- 
vent  has  a  distinctly  lower  temperature  than  other  vents  (203°C),  which  is  even  lower 
than  the  temperature  reported  in  1991  (263°C)  (Von  Damm,  2000). 

3.2.2  Sample  Processing 
3.2.2. 1  Chimney  sulfides 

Hydrothermal  vent  fluids  and  chimney  minerals  were  sampled  at  both  a  high 
temperature  vent  (Bio  9”  vent)  and  a  low  temperature  vent  (K-vcnt).  In  both  cases,  the 
chimney  was  removed  and  hydrothermal  fluids  were  sampled  from  the  remaining  base  of 
the  chimney. 

Hydrothermal  sulfide  samples  were  crushed  between  two  plexiglass  discs  inside  a 
Teflon  bag  using  a  hydraulic  press.  Sulfide  grains  were  collected  between  500  pm  and 
1 .0  mm  sieves  and  mono-mineral  sulfide  phases  were  isolated  by  hand  picking  under  a 
microscope.  For  each  sample,  between  15  and  50  mg  of  mineral  grains  were  picked  in 
order  to  obtain  a  representative  sample.  Care  was  taken  to  select  pure  sulfide  grains 
without  other  mineral  inclusions  but,  in  some  cases,  inclusions  were  unavoidable 
especially  for  sphalerite  (ZnS)  which  includes  other  phases  such  as  pyrite  and 
chalcopyrite.  Representative  splits  of  mineral  separates  were  saved  for  Fe  and  S-isotope 
analysis  and  trace  element  concentrations  (Rouxcl  et  a!.,  in  prep).  Samples  were  weighted 
in  15  mL  Teflon  beakers  and  dissolved  using  5  mL  of  concentrated  HNOj.  After 
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evaporation  on  a  hot  plate  at  60°C,  complete  dissolution  was  achieved  by  a  second 
evaporation  step  using  5mL  of  aqua  regia.  Dry  residue  was  dissolved  in  5  ml  of  6N  HC1 
with  trace  H2O2  by  heating  at  40°C  in  a  closed  vessel. 

When  the  fluids  from  these  sites  were  brought  on  board  the  ship  they  were 
acidified  with  concentrated  HCI  to  pH  ~  I.  The  “dregs”  (insoluble  or  precipitated 
particles  remaining  in  the  Ti-samplers)  arc  recovered  when  the  samplers  are  disassembled 
by  rinsing  with  Milli-Q  water  and  acetone  on  0.45  urn  filters  and  saved  for  chemical 
analysis.  Class- 100  clean  airflow  conditions  were  maintained  during  further  manipulation 
of  all  samples. 

3. 2.2.2  Hydrothermal  fluids 

Vent  fluids  from  the  MIT  collection  were  mixed  with  6N  HCI  to  a  concentration 
of  2  M  Cl'  for  purification  by  anion  exchange  chromatography.  Because  vent  fluids  from 
the  9°N  site  contained  some  precipitate,  12  mL  of  the  homogenized  flu  id -particle  mixture 
was  reacted  overnight  in  a  closed  PFA  capsule  with  3  mL  14N  HNO3  and  3  mL  6N  HCI 
in  order  to  oxidize  and  dissolve  precipitates.  The  solution  was  diluted  to  20  mL  with  a  Cl' 
concentration  of  2  M  for  column  purification. 

3.2. 2.3  Purification  bv  anion  exchange  chromatography 

Samples  were  purified  by  anion  exchange  chromatography  using  a  procedure 
modified  from  that  of  Marechal  et  al.  (Marechal  et  aL,  1999).  Samples  were  loaded  in  2N 
HCI  in  order  to  minimize  the  amount  of  acid  used  and  decrease  the  chance  of  premature 
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elution,  as  the  affinity  of  Zn  for  the  resin  is  highest  at  2N  HCI  (Kraus  and  Moore,  1953). 
Zn  was  eluted  in  0.06N  HCI,  which  may  help  avoid  elution  of  other  elements  (Chapman 
ct  al.,  2006).  Samples  were  evaporated  to  dryness  in  5mL  PFA  capsules,  and  reacted  at 
high  heat  overnight  with  200  pL  concentrated  HNO?  and  100  pL  concentrated  HF  to 
remove  residual  silicates  from  the  sample  and  residual  organics  that  may  have  leached 
from  the  column. 

3.2.3  Isotopic  analysis 

Samples  were  redissolved  in  2%  HNO_t  at  a  Zn  concentration  of  50,  100,  or 
200  ppb  for  isotopic  analysis.  All  samples  were  analyzed  with  100  ppb  Cu  (Ultra 
Scientific,  Lot  #D00204)  for  instrumental  mass  fractionation  correction.  Isotope  ratio 
measurements  were  made  on  an  IsoProbc  multi-collector  ICP-MS  (Thermo  Electron 
Corporation,  formerly  Micromass)  using  an  APEX  Q  inlet  system  (ESI)  with  a  75  pL 
min'1  MicroMist  nebulizer  (Glass  Expansion).  The  optional  desolvating  membrane  was 
not  used  as  we  have  sometimes  observed  large  instrumental  mass  bias  when  using  it. 
Signal  intensity  was  measured  on  masses  60,  63,  64,  66,  67,  and  68  using  Faraday 
collectors  in  fifteen  blocks  of  ten  seconds.  64Nf  was  subtracted  from  the  ,>4Zn  signal  by 
monitoring  the  signal  from  '“Ni*.  Because  we  used  Al  sampling  cones,  the  corrections  for 
Ni  were  insignificant.  All  peaks  were  corrected  for  acid  blanks  using  an  on-peak  zero 
immediately  prior  to  measurement.  All  Zn  isotope  data  is  reported  in  terms  of  deviation 
from  the  “Lyon  JMC”  Zn  standard. 
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Data  was  corrected  for  instrumental  mass  bias  using  the  Cu  internal  spike.  The 
relationship  between  Cu  and  Zn  instrumental  fractionation  has  been  shown  to  differ  from 
predictions  of  the  exponential  mass  bias  law  (Albaredc  and  Beard,  2004;  Marechal  et  al., 
1999).  A  linear  relationship  between  the  natural  log  of  6SCu/*3Cu  and  “Zn/^Zn  in 
standards  was  used  instead  to  determine  the  mass  bias  relationship.  Because  natural 
changes  in  instrumental  mass  fractionation  arc  greatly  reduced  with  the  removal  of  the 
desolvator  from  the  APEX,  modified  standards  containing  H:S04,  Sr,  or  organics  leached 
from  polyelthylcnc  bottles  after  many  months  of  storage  were  sometimes  used  to 
artificially  increase  mass  bias  (Archer  and  Vance,  2004).  All  samples  were  measured  in 
triplicate  and  were  corrected  for  mass  bias  using  Cu.  Sample  standard  bracketing  with 
Cu-corrcetcd  samples  and  standards  was  used  to  obtain  the  final  values, 

Mass  bias  was  also  corrected  on  several  samples  using  a  Zn  double  spike.  A  spike 
mixed  predominately  from  MZn  and  6'Zn  was  used  added  to  samples  for  a  final 
MZn:w,Zn:<wZn:('1'Zn  ratio  of  4.2  :  1.3:  1.6:  1. Isotope  ratios  were  calculated  using  a 
Newton  Raphson  iteration  (Albarede  and  Beard,  2004). 

Many  of  the  MIT  samples  had  a  small  amount  (less  than  a  gram)  of  a  clear 
crystalline  precipitate  in  the  bottom.  These  samples  had  been  filtered  upon  collection,  so 
the  precipitate  must  have  formed  during  sample  storage.  The  composition  of  the 
precipitate  is  unknown,  but  it  could  be  amorphous  silica  or  elemental  sulfur.  This 
precipitate  was  separated  from  lOmL  of  homogenized  fluid  by  centrifugation  and 
rcdissolvcd  in  2%  HNOj  for  elemental  analysis,  Measurements  were  performed  on  a 
PlasmaQuad  2+  quadrupole  ICP-MS  (VG  Elemental,  now  Thermo  Electron)  using  an  In 
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spike  to  correct  for  sensitivity  changes.  In  all  fluids  the  amount  ofZn  in  the  precipitate 
was  less  than  0. 1%  of  the  dissolved  concentration. 

3.3  Results  and  discussion 

3.3.1  Analytical  accuracy  and  precision 

Error  was  estimated  by  calculating  the  pooled  standard  deviation  of  seven 
samples  each  run  on  two  different  days,  combined  with  the  analytical  error  for  five 
process  standards.  The  calculated  error  was  0.069%o  (2o  s.d.).  The  internal  error  averaged 
over  hundreds  of  samples  was  0.067%o  (2 o  s.d.),  indicating  similar  error  for  samples  and 
standards  after  the  chemical  purification  process.  There  was  no  apparent  difference  in  the 
magnitude  of  internal  error  of  triplicate  analysis  between  hundreds  of  samples  run  at  50. 

1 00,  and  200  ppb. 

Data  was  checked  for  polyatomic  or  doubly  charged  interferences  by  comparing 
666Zn  values  to  668Zn  values  (Fig.  1).  With  two  exceptions,  the  data  showed  the  expected 
trend  in  mass  bias  within  error.  These  two  samples  were  run  twice  and  the  discrepancy 
lies  in  the  668Zn  values,  indicating  that  the  66<’Zn  values  are  accurate.  The  values  obtained 
by  correcting  for  instrumental  mass  bias  using  a  double  spike  matched  those  corrected  for 
mass  bias  using  Cu,  indicating  that  both  mass  bias  corrections  are  robust  (Fig.  2). 
Duplicate  samples  taken  with  two  different  pistons  on  the  Ti-sampler  were  analyzed  for 
both  Bio9  and  K-vcnt  at  EPR  9°N  and  their  isotopic  composition  agrees  within  our 
analytical  precision. 
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666Zn 


Figure  3.1.  Comparison  of  6f’*Zn  and  6H’Zn  values  for  all  samples  in  this  study 
demonstrate  the  absence  of  isobaric  interferences.  The  two  outliers  (circled)  were 
measured  twice  with  matching  66flZn  values.  There  appears  to  be  a  small  unexplained 
bias  towards  heavy  b'^Zn  compared  to  the  values  predicted  by  bMZn. 


3.3.2  Bio  9”  high-temperature  vent 

Mineral  samples  at  this  vent  range  from  bw’Zn  =  -0.09%o  to  +0.32%o.  These  values 
represent  real  differences  in  Zn  isotopes,  since  the  range  is  much  larger  than  our 
analytical  uncertainty.  Yet,  there  is  no  clear  relationship  between  bhhZn  and  location  in 
the  chimney  that  would  allow  us  to  determine  what  process  is  responsible  for  the 
observed  differences.  Changes  in  the  fluid  b6f’Zn  values  across  the  chimney  wall,  or 
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Figure  3.2.  Comparison  of  6w’Zn  values  measured  using  either  a  double  spike  method  or 
an  internal  Cu  spike  to  correct  for  instrumental  mass  bias.  Within  the  analytical  error, 
values  provided  by  both  methods  arc  in  agreement. 


isotope  fractionations  during  sulfide  precipitation  (equilibrium  or  kinetic)  could  be 
responsible.  Our  data  does  not  show  the  same  relationship  between  Zn  isotopes  in  co¬ 
existing  sphalerite  and  chalcopyrite  samples  as  observed  in  previous  studies:  chalcopyritc 
666Zn  samples  were  around  0.4%o  lighter  than  nearby  sphalerite  in  samples  from  the 
Alexandria  deposit  in  Russia  (Mason  ct  al.,  2005).  Although  we  find  that  chalcopyrite  is 
lighter  in  three  of  the  four  locations  sampled,  the  maximum  difference  of  Aw'Zn  is  only 
0. 1 8%o.  At  one  location  in  the  chimney  the  sphalerite  was  heavier  than  the  chalcopyritc. 
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Tabic  3.1.  Data  on  all  individual  fluid  and  mineral  samples  measured.  Sample  ID 
is  the  information  written  on  the  sample  bottle,  n  is  the  number  of  triplicate 
analyses  performed  on  different  days.  Error  is  the  internal  standard  deviation  for 
triplicate  analysis  or  the  pooled  standard  deviation  from  two  separate  triplicate 
analyses.  Temperature  and  pH  are  given  as  the  end-member  values  for  the  vent 
area  (Von  Damm,  1984).  Elemental  concentrations  arc  calculated  for  pure 
hydrothermal  fluid  assuming  that  all  the  Mg  in  samples  comes  from  seawater 
entrained  during  sampling,  data  taken  from  Von  Damm,  1984  (EPR  2TN  and 
Guayamas)  from  data  provided  by  Andy  Campbell  (TAG)  and  data  provided  by 
Olivier  Rouxel  (EPR  9°N). 
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This  disparity  may  result  from  differences  in  the  conditions  under  which  the  minerals 
precipitated.  While  the  minerals  studied  by  Mason  et  al.  were  formed  deep  within  the 
hydrothermal  stockwork,  our  samples  were  taken  from  a  chimney  where  minerals  most 
likely  precipitated  much  more  quickly  as  vent  fluids  mixed  with  seawater.  Precipitation 
might  have  occurred  too  quickly  for  isotopic  equilibrium  to  be  established  between  Zn  in 
sphalerite  and  chalcopyrite,  or  might  quantitatively  precipitate  Zn  within  the  chimney 
wall  overwhelming  any  kinetic  or  equilibrium  fractionations  that  did  occur. 

Hydrothermal  11  u ids  sampled  from  the  Bio9”  chimney  have  an  isotopic 
composition  similar  to  the  chimney  minerals.  The  average  Zn  isotope  composition  of  the 
fluids  was  66hZn  =  +0.20%o,  compared  to  +0.  !3%o  for  the  minerals  sampled.  Zn  isotope 
ratios  for  particulate  and  fluid  fractions  in  the  water  samplers  was  the  same  within 
analytical  error.  The  values  are  also  similar  to  existing  measurements  for  the  6w’Zn  of 
basalt.  The  similarity  between  the  bt(,Zn  of  vent  fluids  and  basalt  suggests  that  chemical 
cycling  of  Zn  at  the  site  is  minimal.  The  range  for  all  reported  measurements  of  Zn 
isotopes  in  basalt  is  +0.  19%o  to  +0.30%o  (Archer  and  Vance,  2002;  Chapman  et  al.,  2006; 
Dolgopolova  et  al.,  2006;  Marechal  et  al.,  2000).  Because  of  the  isotopic  similarity 
between  the  basalt  Zn  source  rock  and  the  hydrothermal  fluids,  we  hypothesize  that 
subsurface  precipitation  of  Zn  sulfides  at  this  site  is  minimal. 

3.3.3  K-vent  low-temperature  vent 

There  was  a  large  isotopic  offset  between  the  fluids  (666Zn  =  +1 ,29%o)  and  dregs 
(&66Zn  =  +0.33%o)  in  the  samplers  from  this  site.  Because  the  fluids  were  observed  to  be 
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Figure  3.3A.  Sample  locations  and  6,’llZn  values  of  sphalerite  (regular)  and 
chalcopyrite  (italic)  samples  taken  from  a  chimney  at  the  liio  9”  vent  site  at  9"N 
on  the  East  Pacific  Rise.  (Modified  from  Rouxel  et  al.,  in  prep.) 
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+0.60%o 


Figure  3.3 B.  Sample  locations  and  666Zn  values  for  sphalerite  samples  taken  from 
a  chimney  at  K-vent  at  9°N  on  the  East  Pacific  Rise.  (Modified  from  Rouxel  et  al., 
in  prep.) 
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Figure  3.3C.  6h('Zn  values  of  chimney  sulfides  and  fluids  from  the  high-temperature 
(380°C)  Bio9”  vent  and  the  low-temperature  (203°C)  K  vent  from  the  hydrothermal 
fields  at  9UN  on  the  East  Pacific  Rise.  Mineral  samples  are  aligned  roughly  from  top  to 
bottom. 
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clear  when  they  exited  the  vent,  we  presume  that  the  particulates  formed  in  the  sampler 
rather  than  the  subsurface  fluids.  Assuming  that  the  particles  precipitated  after  sampling, 
the  6f,<’Zn  of  the  original  fluids  are  calculated  by  isotope  mass  balance  as  6(1(’Zn  =  +0.90%o 
and  -t-O.93%0  for  the  duplicate  samples.  The  fractionation  factors  (Aw'Zn)  calculated  from 
the  Rayleigh  equations  calculated  here  are  -0.79%o  and  -0.73%o. 

The  chimney  minerals  at  this  vent  site  range  from  6(,t’Zn  =  +0.57%o  to  +1 . 1 7%o,  a 
greater  range  than  at  the  high-temperature  vent  Bio  9”  and  significantly  heavier  than  the 
predicted  5f>(,Zn  in  the  basalt  source  rock.  Based  on  our  evidence  that  sulfides  are  lighter 
than  the  fluids  from  which  they  precipitate,  we  believe  that  the  presence  of  isotopically 
heavy  chimney  minerals  reflects  their  precipitation  from  isotopically  heavy  vent  fluids. 
The  variation  in  6(,hZn  values  observed  in  these  chimney  minerals  may  result  from 
temporal  variations  in  the  fluid  6hhZn,  or  variations  in  the  precipitation  isotope  effect. 

3.3.4  Worldwide  hydrothermal  fluids 

No  significant  trend  in  fluid  6w’Zn  was  observed  as  a  function  of  concentrations  of 
Zn,  1/Zn,  Cu,  Fe,  H;S,  HS\  fluid  pH,  or  [Fe]/[Mn],  but  a  negative  correlation  with 
temperature  was  observed  (p  <  0.0001)  (Fig.  4).  All  four  of  the  highest  66t'Zn  values 
occurred  in  fluids  with  relatively  low  Cu  concentrations  (s  1  pM)  and  Fc/Mn  ratios  (s 
1.3). 

Several  factors  may  be  responsible  for  differences  in  fluid  6w'Zn,  including 
fractionation  during  phase  separation,  kinetic  or  equilibrium  fractionation  during 
subsurface  precipitation  ofZn  sulfides,  and  subsurface  redissolution  ot  Zn  sulfides. 
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Because  the  total  range  of  fluid  66ftZn  is  an  order  of  magnitude  greater  than  variations  in 
basalt  (Archer  and  Vance,  2002;  Chapman  ct  al.,  2006;  Marcchal  et  al.,  2000),  we  suspect 
that  source  rocks  arc  not  the  source  of  isotopic  variability.  Phase  separation  also  seems 
unlikely  to  fractionate  Zn  isotopes  because  we  do  not  observe  any  relationship  between 
chlorinity  and  6w,Zn. 

We  hypothesize  that  the  subsurface  precipitation  of  light  Zn  sulfides  in  lower- 
temperature  vents  is  the  main  cause  of  isotopic  variation  in  hydrothermal  fluids.  Sulfide 
precipitation  in  the  laboratory  was  accompanied  by  an  isotope  effect  of  A66Zn  =  -0.36%o 
(Archer  et  al.,  2004).  While  the  magnitude  of  this  fractionation  may  vary  under  different 
hydrothermal  conditions,  we  assume  that  sulfides  arc  lighter  than  the  fluid  from  which 
the  precipitate  under  hydrothermal  conditions  as  well.  Temperature  can  have  a  large 
impact  on  sphalerite  solubility  in  hydrothermal  fluids.  Sphalerite  solubility  in  modeled 
hydrothermal  fluids  increases  by  as  much  as  three  orders  of  magnitude  between  200°C 
and  350°C  (Tivey  ct  al.,  1999).  Cooling  may  therefore  be  an  important  factor  in 
governing  the  extent  of  Zn  sulfide  precipitation. 

Fe/Mn  ratios  and  Cu  and  Zn  concentration  data  are  consistent  with  our  hypothesis 
that  Zn  sulfide  precipitation  causes  hydrothermal  fluid  6w>Zn  to  rise.  Fc  is  more  rapidly 
precipitated  when  hydrothermal  fluids  cool  whereas  Mn  is  conservative  during  fluid 
cooling,  so  low  Fe/Mn  ratios  arc  a  good  indicator  of  subsurface  metal  sulfide 
precipitation  (Seewald  and  Seyfricd,  1990).  The  heaviest  6fc6Zn  fluids  also  have  low- 
concentrations  of  Cu  and  Zn.  This  may  be  due  to  differences  in  metal  leaching  from  the 
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source  rock,  or  to  the  precipitation  of  sulfides.  Cu  and  Zn  sulfides  precipitate  from 
hydrothermal  fluids  under  similar  conditions. 

Either  conductive  cooling  or  mixing  with  seawater  could  cause  the  fluid  cooling 
that  we  observe.  The  low  concentrations  of  Mg  in  fluid  samples  from  K-vent  and  Bio9” 
suggest  that  mixing  with  seawater  was  not  a  factor  at  these  sites.  In  fluid  samples  from 
the  MIT  collection.  Mg  in  the  samples  could  come  from  cither  subsurface  mixing  of 
hydrothermal  fluids  with  seawater  or  mixing  during  sampling. 

3.3.5  The  hydrothermal  Zn  isotope  budget 

Zn  in  hydrothermal  vent  fluids  may  impact  on  the  6w>Zn  of  seawater.  The  mid- 
ocean  ridge  hydrothermal  Zn  flux  has  been  estimated  to  be  1300%  of  the  natural  riverine 
flux  (Wheat  ct  al.,  2002).  Much  of  the  Zn  in  hydrothermal  vent  fluids  precipitates  quickly 
as  metal  sulfides  in  the  region  near  the  vent,  but  some  may  remain  in  hydrothermal 
plumes  cither  as  particles  or  dissolved  Zn  (German  ct  al.,  1991).  Estimates  of  ocean 
cycling  time  through  hydrothermal  plumes  arc  several  thousand  years  (Elderfield  and 
Schultz,  1996).  on  the  order  of  the  residence  time  of  the  deep  ocean.  Isotopic  exchange 
between  hydrothermal  and  deep-water  Zn  in  hydrothermal  plumes  could  affect  seawater 
Zn  isotopes.  We  estimate  the  riverine  and  atmospheric  input  of  Zn  to  the  ocean  to  have  a 
b^Zn  isotopic  composition  between  0.09  and  0.29%o.  This  range  encompasses 
measurements  of  Zn  isotopes  in  basalts,  Niger  dust  (Marechal  ct  al..  2000),  an  average  of 
four  reference  ores  (Chapman  ct  al.,  2006).  and  measurements  of  common  anthropogenic 
Zn  products  (S.  John  et.  al..  submitted).  We  see  vent  fluids  as  heavy  as  1.29%o  in  our 
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Figure  3.4.  &M’Zn  as  a  funclion  of  fluid  temperature  and  dissolved  elemental 
concentrations.  Temperature  is  the  highest  recorded  for  that  vent  during  the 
research  cruise,  elemental  concentrations  were  measured  for  each  individual 
sample  (Von  Damm,  1984)  and  extrapolated  to  the  hydrothermal  endmember 
based  on  the  assumption  of  zero  Mg  concentration  in  the  hydrothermal  fluids. 
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study.  When  hydrothermal  fluids  mix  with  seawater,  even  more  light  Zn  may  be  lost  into 
sulfide  precipitates,  so  Zn  in  hydrothermal  plumes  could  be  even  heavier.  We  have 
measured  the  isotopic  composition  of  deep  Pacific  seawater  to  be  Sw’Zn  =  -H).46%o  (John 
et  al„  2005),  heavier  than  our  estimated  riverine  and  atmospheric  input.  Hydrothermal 
processes  may  be  responsible  for  this  difference. 

3.4  Conclusions 

We  present  the  first  data  on  Zn  isotopes  in  modern  hydrothermal  systems  and  the 
first  data  on  Zn  isotopes  in  hydrothermal  vent  fluids.  There  is  a  strong  negative 
correlation  between  the  temperature  of  the  vent  fluids  and  their  b^Zn.  We  hypothesize 
that  this  results  from  the  precipitation  of  isotopically  light  Zn  sulfides  from  iow- 
temperaturc  fluids.  Cu  concentration  data  and  Fe/Mn  data  in  the  fluids  is  consistent  with 
this  hypothesis.  Zn  isotopes  may  be  a  tool  to  help  us  understand  the  plumbing  and 
chemistry  of  hydrothermal  sy  stems.  Hydrothermal  systems  may  also  be  an  important 
factor  in  setting  the  b66 Zn  of  the  deep  ocean. 
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Chapter  4 


Zinc  isotope  fractionation 
associated  with  two  separate 
uptake  pathways  in  a  marine 
diatom 


When  micronutricnts  such  as  zinc  are  transported  across  a  cell  membrane,  some 
isotopes  may  be  taken  up  more  quickly  than  others.  This  biological  fractionation  is 
recorded  as  deviations  in  the  ratio  of  66Zn/64Zn  in  environmental  samples  such  as 
seawater  (Bcrmin  et  al.,  2006),  marine  sediments  (Pichat  et  al.,  2003),  and  marine  life 
(Marechal  et  al.,  2000).  Measuring  the  magnitude  of  zinc  isotope  fractionation  by  marine 
phtoplankton  is  crucial  to  interpreting  these  natural  signals.  Here  we  distinguish  between 
cellular  Zn  acquired  through  separate  high-affinity  and  low-affinity  Zn  transport  systems. 
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and  calculate  the  isotopic  fractionation  associated  with  each  pathway.  We  have  calculated 
the  biological  isotope  effect  for  high-  and  low-affinity  transport  to  be  A66Zn  =  -0.2%o  and 
-0.8%o,  respectively,  in  the  ecologically  important  marine  diatom  Thalassiosira  oceunica. 
We  have  discovered  these  transport- related  isotope  effects  by  growing  cultures  at  a  range 
of  Zn  concentrations,  and  distinguishing  between  isotopically  light  intracellular  Zn  and 
the  isotopically  heavy  Zn  which  adsorbs  to  the  ceil  exterior.  Consequently,  this  study  is 
the  first  to  describe  a  molecular  basis  for  metal  isotope  fractionation  during  transport.  A 
switch  in  the  predominance  of  high-affinity  and  low-affinity  transport  occurs  in  many 
species  of  phytoplankton  over  an  oceanographically  relevant  range  of  Zn  concentrations 
(Sunda  and  Huntsman,  1992;  Sunda  and  Huntsman,  1998),  suggesting  that  the  isotope 
fractionations  seen  here  may  be  common  in  the  ocean.  By  demonstrating  that 
fractionation  values  can  be  tied  to  specific  transporters,  we  show  that  natural  variations  in 
transition  metal  isotope  ratios  can  be  interpreted  in  terms  of  the  underlying  molecular 
biological  processes  from  which  they  arise.  By  studying  the  transport-associated 
fractionation  of  other  metals  and  in  other  organisms,  transition  metal  isotopes  have  the 
potential  to  address  many  quest  ions  of  oceanographic,  environmental,  and  medical 
importance. 
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4.1  Introduction 


Zinc  is  an  essential  biological  nutrient  in  the  oceans  where  the  concentration 
distribution  of  Zn  is  controlled  largely  by  phytoplankton  Zn  uptake  and  remineralization 
(Bruland  and  Lohan,  2003).  Biological  Zn  isotope  effects  of  up  to  0.5%o  in  b^Zn  have 
been  observed  in  land  plants  (Weiss  et  al.,  2005)  and  phytoplankton  (Gelabert  et  al., 
2006)  and  are  substantial  compared  to  a  total  reported  range  for  natural  6w’Zn  values  of 
about  3.5%o  (Wilkinson  et  al.,  2005),  where: 
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The  biological  fractionation  ofZn  isotopes  has  been  invoked  to  interpret  natural  Zn 
isotope  signals  in  seawater  (Bcrmin  et  al.,  2006),  seafloor  sediments  (Pichat  et  al.,  2003), 
marine  particles,  and  manganese  nodules  (Marechal  et  al.,  2000).  By  accurately 
measuring  the  magnitude  of  biological  Zn  isotope  fractionation,  changes  in  the  isotopic 
composition  of  Zn  in  seawater  or  ancient  sediments  can  be  directly  related  to  the  extent 
of  biological  Zn  uptake  in  surface  waters.  The  stable  isotopes  of  C  and  N  have  been  used 
for  many  years  to  study  environmental  and  biological  processes;  Zn  isotopes  may  have  a 
similar  potential  to  help  us  understand  present  and  ancient  oceans. 
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4.1.1  The  kinetics  of  high-affinity  and  low-affinity  transport 

We  have  grown  the  marine  diatom  Thalassiosira  oceanic  a  in  EDT A- buffered 
seawater  media  at  free-Zn  concentrations  spanning  the  natural  range  of  the  marine 
environment,  from  10'1*  M  in  the  nutrient-poor  central  ocean  gyres  (Bruland.  1989)  up  to 
10"  M  in  anthropogenically  impacted  coastal  waters  such  as  Narragansctt  Bay  ( Kozclka 
and  Bruland.  1998).  In  order  to  distinguish  between  extracell ularly  adsorbed  Zn  and  Zn 
that  has  been  transported  to  the  cell  interior,  filtered  cells  from  each  culture  were  split 
into  two  fractions  and  washed  either  with  Zn-ffee  seawater  or  with  an  EDTA-oxalatc 
wash  that  removes  adsorbed  metals  (Tang  and  Morel,  2006;  Tovar-Sanchez  et  al.,  2003; 
Tovar-Sanchez  et  al.,  2004).  By  removing  extracellular  metals,  we  can  isolate  the 
fractionation  ofZn  associated  with  transport  across  the  cell  membrane.  The  internalized 
specific  Zn  uptake  rate  measured  in  our  experiments  shows  a  sigmoidal  relationship  to  Zn 
concentration,  consistent  with  the  shift  between  low-affmity  and  high  affinity  zinc 
transporters  that  has  been  observed  to  occur  in  marine  diatoms  (Sunda  and  Huntsman, 
1992)  (Fig.  2).  This  sigmoidal  uptake  curve  has  been  described  by  the  cumulative  uptake 
of  high  and  low  affinity  transporters  governed  by  the  equation: 


V  = 


JZ«2+] 


\Z>rl+KM 


+  AlZn2*}. 


(2) 


The  uptake  rate  ( V)  is  the  cumulative  effect  of  two  terms  representing  the  high-affinity 
and  low-affinity  uptake  patways.  High-affinity  Zn  uptake  predominates  at  low  [Zn2  ]  and 
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is  governed  by  the  first  term,  following  classical  Michaclis-Menten  kinetics  dependant  on 
the  maximal  uptake  rate  (f'max)  and  the  Michaelis-Menten  constant  (Am)-  At  high  Zn 
concentrations,  low-affinity  transport  dominates  with  an  uptake  rate  related  to  [Zn‘  ] 
concentration  by  a  constant  (A).  The  decreasing  contribution  of  high-affinity  uptake  at 
higher  [Zn2+]  can  result  either  from  saturation  or  down-regulation  of  the  high-affinity 
transporters. 

These  high-affinity  and  low-affinity  transport  systems  for  Zn  have  been  identified 
in  many  organisms,  including  fish  (Glover  et  al.,  2003;  Qiu  and  Hogstrand,  2005),  yeast 
(Zhao  and  Eide,  1996),  and  many  phytoplankton  (Sunda  and  Huntsman,  1992;  Sunda  and 
Huntsman,  1998).  ZIP-family  proteins,  many  of  which  are  specific  Zn  transporters,  have 
been  found  in  bacteria,  archaca,  and  many  eukaryotes,  and  putative  ZIP  homologues  have 
been  identified  in  the  genome  of  the  diatom  T.  pseudonana,  Several  other  eukaryotic 
phytoplankton  have  similar  Zn  transporter  affinities  and  exhibit  a  similar  switch  between 
high-affinity  and  low-affinity  transport  systems,  including  three  species  of  diatoms,  the 
coccolithophore  E.  huxleyii ,  and  the  alga  Chlamydomonas  (Sunda  and  Huntsman,  1992; 
Sunda  and  Huntsman,  1 998).  Similar  Zn  transport  kinetics  suggests  analogous  Zn  isotope 
fractionations  by  these  different  species. 

4.2  Methods 
4.2.1  Culture  growth 

The  diatom  T.  oceanica  (CCMP  1005)  was  maintained  in  acid-cleaned  polycarbonate 
bottles  at  20°C  under  85  p Einsteins  m'2  s'1  constant  light.  Media  was  prepared  from 
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Figure  4. 1 .  The  bioiogicai  isotope  effect  during  Zn  uptake  is  different  for  high-affinity 
and  low- affinity  transport.  Suggested  Zn  transport  by  high -affinity  (dark  grey)  and  low- 
affinity  (light  grey)  transporters  under  low-Zn  (A)  and  high-Zn  (B)  conditions.  When  Zn 
concentrations  are  low  (c.g.  central  oceans),  high-affinity  transporters  arc  uprcgulated 
and  account  for  most  Zn  transport  with  an  isotope  effect  of  A66Zn  =  -0.2%<>.  Under  high- 
Zn  concentrations,  Zn  uptake  occurs  primarily  through  a  low-affinity  transport  pathway 
with  an  isotope  effect  of  A^Zn  =  -0.8%o. 


filtered  Sargasso  surface  seawater  collected  with  a  towed  trace-element  clean  Fish 
modified  with  a  rigid  rubber  tube  in  place  of  the  PVC  pipe  used  for  support  (Vink  ct  al., 
2000)  amended  with  major  and  trace  element  nutrients  and  vitamins  (Sun da  and 
Huntsman,  1995)  with  10  pM  Fe,  25  nM  Co,  and  varying  amounts  of  Zn.  Added  Zn  was 
a  laboratory  standard  prepared  with  Zn  from  the  Canadian  Electrolytic  Zinc  (CEZ) 
corporation.  This  Zn  has  a  similar  isotopic  composition  to  natural  Zn  (S.  John,  submitted 
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to  Chemical  Geology)  so  small  amounts  ofZn  contamination  would  not  have  a  large 
effect  on  666Zn.  Seawater  was  microwave  sterilized  and  all  nutrients,  vitamins,  and  trace- 
elements  were  sterilized  by  syringe  filtration  with  an  acid-cleaned  rubber/silicone- free 
syringe  (NormJect,  HSW)  through  an  acid-cleaned  0.4  pm  polycarbonate  filter 
(Nuclepore)  and  added  to  the  seawater  after  it  had  cooled.  0. 1%  inoculums  by  volume 
for  experiments  were  taken  from  cultures  in  log-phase  growth  that  were  maintained  in 
media  with  a  Zn  concentrations  within  one  log  unit  of  experimental  conditions.  Growth  in 
1L  polycarbonate  bottles  at  maintenance  conditions  was  monitored  by  measuring  in-vivo 
Chl«  fluorescence  with  a  Synergy  HT  Microplate  reader  (BioTck).  Specific  Zn  uptake 
rate  was  calculated  by  multiplying  the  log- phase  growth  rate  by  the  Zn:P  ratio  of  the  final 
cultures. 

4.2.2  Processing  cultures  for  isotopic  analysis 

Cultures  were  harvested  near  the  end  of  log-phase  growth  and  filtered  onto  a 
1  pm  polycarbonate  filter  (Whatman  Nuclepore).  Each  culture  was  split  into  two  equal 
parts,  one  of  which  was  rinsed  once  with  mctal-frcc  seawater  and  the  other  was  rinsed 
with  an  oxalatc-EDTA  reagent  to  remove  adsorbed  extracellular  metals  (Tovar-Sanchcz 
ct  al„  2003;  Tovar-Sanchcz  et  al„  2004)  modified  to  pH=7  for  two  5-minute  rinses  to 
ensure  complete  removal  of  extracellular  ferric  oxides  (Tang  and  Morel,  2006).  Cells  and 
filters  were  transferred  to  -5  mL  quartz  beakers  and  organics  were  partially  oxidized  be 
reacting  for  several  hours  with  2  mL  warm  HNOj  and  200  pL  HiO:.  When  the  liquid  had 
evaporated  the  remaining  organics  were  oxidized  by  combusting  samples  for  8  hours  at 
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450°C.  Combusted  material  was  dissolved  in  warm  6N  HC1  for  column  purification  of  Fc 
and  Zn,  small  subsamplcs  were  diluted  for  measurement  of  Zn:P  ratios. 

Zn  concentrations  were  measured  by  isotope  dilution  ICP-MS  with  a  6KZn  spike.  P 
concentrations  were  performed  colorimctrically  according  to  standard  procedures  for 
measurements  of  P  in  natural  seawater  (Strickland  and  Parsons,  1972).  Zn  and  Fc  were 
purified  by  anion  exchange  chromatography  (Marechal  ct  al„  1999)  using  small  column 
volumes  to  reduce  metal  contamination  from  acids  (Archer  and  Vance,  2004)  and  the  Zn 
fraction  was  eluted  with  0.1%  HCI  to  avoid  co-clution  of  other  elements  (Chapman  et  al„ 
2006).  Samples  were  evaporated,  then  reacted  overnight  in  5  mL  PFA  capsules  (Savillex) 
at  high  temperature  with  200  pL  HF  and  100  pL  HjCF  to  drive  off  silicates  remaining 
from  the  diatom  frustules  and  organics  that  may  have  leached  off  the  column.  Samples 
were  evaporated  to  dryness,  evaporated  again  with  100  pL  I6N  HNOi  to  drive  off 
remaining  F'.  and  dissolved  in  2%  HNO.i  for  isotopic  analysis. 

4.2.3  Isotope  analysis 

Samples  were  analyzed  for  Zn  isotopes  on  an  IsoProbe  multi-collector  ICP-MS 
(GV  Instruments,  formerly  MicroMass)  equipped  with  an  Apex  Q  inlet  system  without 
the  optional  dcsolvator.  Samples  were  prepared  with  50  ppb  Cu  and  either  25,  50,  or 
100  ppb  Zn.  Signal  was  monitored  on  masses  60,  63,  64,  65,  66,  67,  and  68.  Corrections 
for  instrumental  mass  bias  were  made  by  monitoring  the  65Cu/6,Cu  ratio  of  the  Cu  spike 
in  samples  and  using  interspersed  standards  to  establishing  a  linear  mass  bias  relationship 
between  the  natural  log  65Cu/63Cu  and  the  natural  log  of  66Zn/’4Zn  (Marechal  ct  al„ 
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1999).  MNi'  was  subtracted  from  wZn’  by  monitoring  60Ni'  and  applying  the  Cu  mass 
bias  correction,  although  these  corrections  were  insignificant.  6™Zn  values  showed  twice 
the  fractionation  compared  to  6w’Zn  values  and  indicate  the  absence  of  polyatomic 
interferences.  Extensive  error  analysis  based  on  hundreds  of  samples  was  performed  in 
conjunction  with  a  study  of  Zn  isotopes  in  hydrothermal  systems  (S.  John,  in 
preparation).  External  reproducibility  is  0.07%o  (2a  s.d.)  independent  of  the  Zn 
concentration  in  the  analysis. 

4,3  Results  and  discussion 
4.3.1  Zn  uptake  kinetics 

We  have  modeled  our  specific  Zn  uptake  rate  data  as  the  sum  of  high-affinity  and 
low-affinity  uptake,  assuming  that  nearly  all  of  the  uptake  at  10'*'"  M  Zn*'  occurs  through 
low-affinity  transport  (Fig.  2).  According  to  Eqn.  2,  Zn  uptake  through  the  low-affinity 
transport  system  decreases  linearly  with  Zn  concentration.  Any  uptake  in  excess  of  the 
predicted  low -affinity  Zn  uptake  is  assumed  to  have  occurred  through  the  high-affinity 
Zn  transport  pathway.  Accordingly,  nearly  100%  of  Zn  uptake  occurs  though  the  low- 
affinity  transport  system  at  1 0‘4  M  Zn‘  and  above,  there  is  a  transition  between  the 
predominance  of  these  two  uptake  path  ways  between  10’11>  "  and  10’)5  M  Zn'  ,  and  low- 
affinity  transport  accounts  for  less  than  1 0%  of  total  Zn  uptake  when  the  concentration  of 
Zn:  is  10'"  M  and  below  (Fig.  3). 
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4.3.2  Biological  Z n  isotope  fractionation 

The  Zn  isotope  ratio  in  our  diatoms  is  correlated  to  the  switch  in  Zn  transport 
systems  (Fig  3).  The  values  of  Aw’Zn  (6fl<iZn<jiaioms  -  &f’hZnm«jia)  changes  rapidly  over  the 
same  range  of  Zn  concentrations  where  the  Zn  uptake  switches  from  predominantly  high- 
affinity  uptake  to  predominantly  low-affinity  uptake  when  Zn:  is  between  t  O' 10  5  M  and 
l  O'1’ >  M.  Aw’Zn  is  around  -0.2%o  at  the  lowest  Zn  concentrations  and  around  -0.8%o  at  the 
highest  Zn  concentrations,  with  intermediate  values  when  there  is  significant  uptake 
through  both  the  high-affinity  and  low-affinity  transport  systems.  We  have  averaged 
values  of  A!,,’Zn  at  the  three  lowest  Zn  concentrations  and  the  two  highest  Zn 
concentrations  resulting  in  the  first  estimates  of  the  isotopic  fractionation  factors  for  Zn 
transport,  with  Aw’Zn  -  -0.2%o  for  high-affinity  transport  and  Aw'Zn  =  -0,8%o  for  low 
affinity  transport. 

We  hypothesize  that  the  fractionation  observed  during  high-affinity  transport 
probably  does  not  occur  at  the  transporter  itself,  but  occurs  because  the  6(,<’Zn  of  frec-Zn 
arriving  at  the  transporter  is  lighter  than  Zn  in  the  bulk  media.  At  equilibrium,  free-Zn  is 
-0. 16%o  lighter  than  EDTA-bound  Zn  in  aqueous  solution  (Ban  ct  al.,  2002).  An 
additional  isotope  effect  will  be  conferred  by  the  difference  in  diffusivity  between  the 
different  Zn  isotopes.  MZn’'  will  diffuse  towards  the  cell  more  quickly  than  66Zn2',  and 
the  difference  in  aqueous  diffusivity  of  free  Zn  is  equivalent  to  a  preference  for  light  Zn 
uptake  by  A'^Zn  =  -0.06%o  (Rodushkin  et  al.,  2004).  Together,  the  free  Zn/Zn-EDTA 
equilibrium  isotope  effect  and  the  difference  in  isotope  diffusivity  can  account  for  the 
entire  -0.2%o  fractionation  observed  under  high-affinity  transport.  Zn  concentrations  in 
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Figure  4.2.  Zn  uptake  rates  show  a  concentration-dependant  switch  in  the  predominance 
of  high-affinity  and  low-affinity  Zn  transport.  Specific  Zn  uptake  rates  measured  for 
washed  (O)  T.  oceanica  cells  at  different  concentrations  of  Zn2',  and  the  predicted 
specific  Zn  uptake  rate  by  the  low-affinity  transport  system  alone  ( — ). 
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the  open  ocean  are  nearly  low  enough  to  limit  phytoplankton  growth  (Ell wood  and  Van 
den  Berg,  2000;  Lohan  et  al.,  2002),  so  it  stands  to  reason  that  high-affinity  transporters 
efficiently  acquire  all  Zn  diffusing  to  the  transporters  regardless  of  isotope.  By  contrast, 
the  isotope  effect  associated  with  the  low-affinity  uptake  pathway  is  much  larger.  This 
isotope  effect  could  result  from  preferential  retention  of  the  heavier  Zn  isotopes  to  low- 
affinity  transport  binding  site,  or  a  more  rapid  translocation  of  light  Zn  isotopes  across 
the  cell  membrane  such  as  would  result  from  faster  diffusion  through  an  ion  channel. 

When  cells  are  not  washed  to  remove  extracellular  metals,  external  Zn  can  have  a 
major  impact  on  the  total  cellular  Aw’Zn.  In  this  experiment,  the  A66Zn  of  unwashed  cells 
ranged  from  -0.05%«  at  the  lowest  Zn  concentration  to  +0.38%<>  at  the  highest  Zn 
concentration.  The  contribution  of  internal  Zn  isotopes  can  be  subtracted  to  calculate 
values  of  A""Zn  for  the  external  Zn  pool  alone.  We  have  found  that  external  Aw’Zn 
increases  almost  linearly  from  +0.09%i>  at  10'"  5  Zn2'  to  +0.52%<>  at  10'8  SM  Zn24.  These 
results  are  consistent  with  other  studies  showing  that  heavy  Zn  isotopes  arc  preferentially 
adsorbed  or  precipitated  onto  diatom  exteriors  (Gelabert  ct  al.,  2006).  In  previous 
experiments  (unpublished  data)  we  have  found  the  amount  of  externally  bound  Zn  and 
the  consequent  values  of  cell  total  Aw>Zn  arc  be  highly  variable  between  different 
experiments.  Cell-surface  precipitation  may  be  a  major  contribution  to  total  Zn  when 
cells  are  grown  in  EDTA-buffered  seawater  media  because  of  high  concentrations  of  total 
Zn  and  high  rates  of  Fe-oxide  precipitation  (Tang  and  Morel,  2006).  In  fact,  simply 
reducing  the  concentration  of  free  Fe  in  the  media  can  reduce  the  levels  of  extracellular 
Zn  to  a  fraction  of  the  intracellular  concentration  (Tang  and  Morel,  2006).  Laboratory 
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Figure  4.3.  The  transition  between  high-affinity  and  low-affinity  transport  is 
conincident  with  a  change  in  the  biological  isotope  effect.  Predicted  specific  Zn 
uptake  through  the  low-affinity  transport  system  as  a  percentage  of  the  total 
specific  Zn  uptake  rate  {-  -  -)  and  the  offset  between  &MlZn  of  T.  ocean ica  cells 
and  the  b66Zn  of  the  seawater  media  they  were  grown  in  (A(,flZn)  at  different 
concentrations  of  free  Zn2+  (O).  The  shift  in  A(,f’Zn  occurs  within  a  similar  range 
of  Zn  concentrations  to  the  shift  between  the  predominance  of  high-affinity  and 
low-affinity  uptake  indicating  that  each  transport  pathway  has  a  unique  isotope 
effect.  Error  bars  represent  external  reproducibility  (2a  s.d.). 
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Figure  4.4.  Unwashed  cells  are  isotopically  heavier  than  washed  cells.  A  comparison 
cells  washed  to  remove  extracellular  precipitates  (O)  with  unwashed  cells  { • )  shows 
that  extracellular  precipitates  are  isotopically  heavy  and  can  greatly  influence  the  total 
cell  isotope  ratio  if  not  removed. 
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studies  are  potentially  fraught  with  errors  associated  with  extracellular  precipitation,  a 
phenomenon  that  is  not  likely  to  be  important  in  all  but  the  most  polluted  or  concentrated 
marine  environments. 

4.4  Conclusions 

Our  experiments  distinguish  between  several  distinct  processes  that  may 
contribute  to  the  total  observed  signal  for  “biological  fractionation”  in  the  oceans.  Cell- 
surface  precipitation  may  not  always  be  a  significant  contribution  to  cellular  Zn,  but 
when  precipitation  docs  occur  we  expect  a  preference  for  heavy  Zn  to  adsorb  onto  the 
cell  surfaces.  An  isotopic  fractionation  between  free  Zn  and  Zn  bound  to  organic  ligands, 
as  we  have  observed  with  Zn  binding  to  EDTA,  may  also  be  important  in  the  oceans 
where  most  of  the  dissolved  Zn  is  organically  complexed  (Bruland,  1989;  Ell  wood  and 
Van  den  Berg.  2000).  Finally,  we  expect  the  isotopic  composition  of  phytoplankton’s 
internal  Zn  pool  to  be  set  by  the  isotope  effects  of  high-affinity  and/or  low-affinity  uptake 
depending  on  the  ambient  concentrations  of  free  Zn. 

As  more  data  becomes  available  on  Zn  isotopes  in  the  oceans,  understanding  how 
phytoplankton  fractionate  Zn  isotopes  will  be  crucial  to  interpreting  natural  signals. 
Dissolved  666Zn  variations  of  0.2%o  have  been  measured  between  surface  and  deeper  NE 
Pacific  seawaters  (Bcrmin  ct  al.,  2006),  attributed  to  a  biological  drawdown  of  Zn.  With  a 
knowledge  of  the  molecular  mechanisms  responsible  for  isotopic  fractionation.  Zn 
isotopes  may  eventually  be  used  to  trace  nutrient  supply,  dust  deposition,  Zn/N 
utilization,  metal  spcciation,  and  other  important  biogcochemical  processes  in  the  oceans. 
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Chapter  5 


Zinc  isotopes  in  seawater  and 
natural  marine  plankton 


The  zinc  distribution  in  the  oceans  is  largely  controlled  by  biological  uptake  and 
remineralization.  Phytoplankton  in  culture  fractionate  Zn  isotopes,  so  measurements  of 
Zn  isotope  ratios  may  be  a  valuable  tool  for  tracing  biological  processes  in  the  oceans. 
Here,  I  report  data  on  Zn  isotopes  in  seawater  and  natural  phytoplankton  samples  in  order 
to  better  understand  the  marine  distribution  of  Zn  isotopes  and  to  better  understand  how 
natural  biological  processes  fractionate  Zn  isotopes.  A  new  isotope  dilution  method  tor 
the  measurement  of  Zn  isotope  concentrations  in  seawater  has  been  developed  that  is 
compatible  with  high-precision  isotope  ratio  measurements.  Methods  have  also  been 
developed  for  the  separation  and  purification  of  Zn  from  large  volumes  of  seawater  tor 
isotopic  analysis.  Zn  isotope  ratios  were  measured  in  samples  trom  five  depths  in  the 
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north  Pacific,  and  in  deep  waters  from  the  north  Atlantic  (0.4  -  0.5  %o).  There  appears  to 
be  a  trend  towards  lighter  values  of  6w>Zn  in  surface  waters  that  correlates  well  with  Zn 
concentrations.  Based  on  this  isotopic  evidence,  we  suggest  that  much  of  the  Zn 
transported  to  the  deep  ocean  is  not  isotopically  light  “primary”  biological  Zn,  but  instead 
results  from  the  adsorption  of  isotopically  heavy  Zn2  onto  falling  biological  material. 
Phytoplankton  samples  from  the  Peru  Upwclling  Region,  as  well  as  several  other 
locations  across  the  globe,  show  that  most  natural  phytoplankton  have  6ft6Zn  between  0%o 
and  +0.6%o.  There  is  no  significant  change  in  the  6w‘Zn  of  natural  plankton  over  large 
changes  in  nutrient  concentrations  in  the  Peru  Upwelling  Region.  This  suggests  cither 
that  the  isotope  effect  for  Zn  uptake  in  this  region  is  very  small,  or  that  nutrient  supply  is 
predominantly  vertical,  rather  than  horizontal. 

5.1  Introduction 

5.1.1  Zinc  concentration  distribution  in  the  oceans 

Zn  has  a  “nutrient- like”  vertical  concentration  profile  in  the  oceans  (Fig.  1).  The 
shape  of  this  profile  is  similar  to  Si,  with  a  deeper  regeneration  maximum  than  N  and  P. 
As  with  other  nutrients,  the  concentration  of  Zn  increases  from  the  Atlantic  to  the 
Southern  Ocean  to  the  Pacific.  The  Zn  maximum  increases  from  around  2  nM  in  the 
North  Atlantic  to  6  nM  in  the  Southern  Ocean  to  9  nM  in  the  North  Pacific  (Bruland  ct 
al.,  1978;  Coalc  et  al.,  2005;  Ellwood  and  Van  den  Berg,  2000;  Fitzwatcr  ct  al.,  2000).  In 
the  surface  ocean,  concentrations  of  Zn  can  be  depleted  by  several  orders  of  magnitude 
with  reported  concentrations  as  low  as  0.04  nM  (Lohan  et  al.,  2002)  in  the  surface  North 
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Figure  5,1.  The  vertical  concentration  profiles  of  Zn  ( - ),  Si  (-  -  -),  and  N  (-  -)  in  the 

North  Pacific  Ocean.  Zn  data  from  {Bruland,  1980),  N  and  Si  data  from  GEOSECS 
(GEOSECS  Operations  Group,  1987). 
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Pacific  and  0.06  nM  in  the  surface  North  Atlantic  (Bruland  and  Franks,  1983). 
Throughout  much  of  the  oceans,  Zn  is  mostly  present  complexed  to  organic  ligands  rather 
than  as  inorganic  “free -Zn”  (Bruland  and  Lohan,  2003). 

5.1,2  Zinc  isotope  analysis  in  seawater 

Zinc  isotopes  have  been  measured  in  seawater  by  several  different  methods.  Zn 
has  been  extracted  from  seawater  by  adsorption  onto  Chelex-100  (Bennin  ct  at.,  2004a; 
Bcrmin  ct  al.,  2004b)  and  co-precipitation  with  Mg(OH):  (Bcrmin  ct  ah,  2006;  John  ct 
ah,  2005).  The  possible  effects  of  isotope  fractionation  during  Zn  extraction  can  be 
compensated  for  by  quantitative  recovery  of  the  Zn,  or  by  adding  a  Zn  double  spike  to  the 
sample  before  concentration.  Instrumental  mass  bias  has  been  evaluated  by  cither  a  Zn 
double  spike  or  an  internal  Cu  spike.  Different  methods  of  isotope  analysis  have  different 
strengths.  A  Zn  double-spike  method  does  not  require  that  Zn  recovery  be  quantitative 
and  allows  for  the  mass  bias  correction  to  be  done  using  Zn,  rather  than  Cu  which  may 
have  a  slightly  different  instrumental  mass  bias.  By  avoiding  a  double  spike,  it  is  possible 
to  compare  “ZnA’Zn  values  with  <’sZn/64Zn  in  order  to  evaluate  the  presense  of  iso  baric 
interferences. 

Using  Chelex  extraction,  the  S6f’Zn  of  English  Channel  water  was  measured  to  be 
+0,46%o  using  a  Cu-correction  and  +0.38%o  using  a  double-spike  approach  (Bcrmin  ct  ah, 
2004a).  Mg(OH)2  co- precipitation  and  a  Cu  correction  was  used  to  measure  a  bM'Zn  = 
+0.46%o  for  deep  North  Pacific  waters  (John  ct  ah,  2005).  A  profile  of  Zn  isotopes  from 
the  North  Pacific  using  a  Zn  double-spike  method  includes  lighter  66f’Zn  values  ranging 
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from  -0.  t%o  to  +0,2%o.  This  study  expands  the  range  of  seawater  samples  in  which  Zn 
isotopes  have  been  measured  and,  by  employing  a  different  analytical  methodology, 
provides  an  important  check  on  previous  measurements. 

9.1,3  Zinc  isotopes  in  marine  plankton 

Zinc  isotopes  have  previously  been  measured  in  biological  material  such  as  plants 
(Weiss  et  al.,  2005),  tissue  from  several  marine  species  (Marechal  et  al.,  2000),  and 
cultured  phytoplankton  (Gelabert  et  al.,  2006).  The  isotope  effect  of  extracellular  Zn 
precipitation  and  biological  transport  has  been  determined  for  diatoms  grown  in  culture 
(Chapter  4).  Diatoms  preferentially  acquired  light  Zn  isotopes  from  the  surrounding 
media,  and  the  magnitude  of  this  fractionation  depended  on  the  cellular  machinery  used 
to  transport  the  Zn.  The  isotope  effect  for  transport  though  a  high-affinity  pathway  was 
measured  to  be  A^Zn  =  -0.2%o  and  for  low-affinity  uptake  A(,<’Zn  -  -0.8%o.  At  the  low  Zn 
concentrations  that  predominate  in  the  open  oceans,  high-affinity  transport  will  dominate 
total  Zn  uptake.  In  coastal  regions  where  Zn  concentrations  are  much  higher,  low-afTmity 
transport  may  contribute.  The  Zn  isotope  ratio  for  natural  phytoplankton  will  be  equal  to 
the  S66Zn  of  the  surrounding  seawater  plus  the  biological  isotope  effect.  By  measuring  Zn 
isotopes  in  natural  phytoplankton  and  seawater,  we  can  better  understand  how  Zn 
isotopes  may  be  used  as  a  tracer  for  biological  processes  in  the  modern  and  ancient 
oceans. 
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5.2  Methods 


5.2.1  Trace-element  clean  laboratory  practices 

Class- 100  clean  airflow  conditions  were  maintained  during  manipulation  of  all 
samples  subsequent  to  collection.  Water  was  distilled  in  a  borosilicatc  glass  still  (Coming 
MegaPure)  to  remove  metal  contaminates,  and  all  acids  were  triply  or  quadrupally 
distilled  in  Vycor®  glass  to  remove  contaminants. 

5.2.2  Sample  collection 

Seawater  samples  from  the  North  Pacific  (47°N,  170.5°W)  were  collected  aboard 
the  R/V  Kilo  Moana  in  2003.  Samples  from  the  North  Atlantic  ( I6.5°N,  59°W)  were 
collected  aboard  the  R/V  Endeavor  in  2005,  and  samples  from  the  Peru  Upwelling 
Region  were  collected  aboard  the  R/V  Knorr  in  2006. 

Seawater  samples  from  the  North  Pacific  were  collected  using  Teflon-coated  GO- 
FLO  samplers  (General  Oceanics)  attached  to  a  Kevlar  wire  for  all  samples  at  1000  m 
and  above.  GO-FLO  samplers  were  attached  to  a  stainless  steel  CTD  cable  to  collect 
waters  below  1000  m.  Seawater  samples  from  the  North  Atlantic  were  collected  using 
individual  MITESS  samplers  modified  to  hold  1L  bottles  and  attached  to  a  hydrowirc 
with  a  VANE  to  keep  the  sampler  upstream  of  the  wire  (Bell  ct  al„  2002;  Bergquist, 
2004).  Seawater  samples  from  the  Peru  Upwelling  Region  were  collected  using  either  a 
MITESS/VANE  sampler  or  Teflon-coated  Niskin  bottles  attached  to  a  plastic-covered 
frame  and  suspended  from  a  Kevlar  hydrowire.  All  seawater  samples  were  filtered 
through  a  0. 2pm  or  0. 4pm  polycarbonate  filter  and  acidified  to  pH  2  immediately  after 
collection  with  HC1. 
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Plankton  tows  were  collected  using  a  50  pm  mesh  nylon  net  with  no  metal  fittings 
attached  to  a  Kevlar  hydrowirc  with  plastic  rope  (Sea  Gear  Corp.).  Tows  were  conducted 
in  the  upper  mixed  layer  above  the  chlorophyll  maximum.  Tows  ranged  in  time  from  five 
minutes  to  thirty  minutes  depending  on  the  productivity  of  the  region.  Concentrated 
plankton  from  the  cod  end  was  mixed  vigorously  and  split  into  to  two  equal  portions 
before  filtration.  Plankton  were  filtered  onto  1  pm  polycarbonate  filters  and  washed 
either  with  clean  surface  seawater  or  a  plankton  wash  designed  to  remove  extracellular 
metals  (Tang  and  Morel,  2006;  Tovar-Sanchez  ct  al.,  2003;  Tovar-Sanchez  et  ah,  2004) 
followed  by  a  seawater  rinse,  filtered  phytoplankton  were  frozen  for  storage.  When 
cuphausiidac  (krill)  were  found  in  filtered  plankton  tows,  they  were  picked  out  with 
Teflon  tweezers  and  washed  to  remove  extracellular  metals  in  the  same  manner  described 
for  bulk  plankton. 

5.2.3  Zn  concentration  in  seawater  measurements 

Zn  concentration  in  seawater  was  determined  by  isotope  dilution  I  CP- MS  after 
coprecipitation  with  Mg(OH)2  (Wu  and  Boyle,  1998)  and  purification  by  anion  exchange 
chromatography.  A  known  weight  of  each  sample  (approximately  1 5  mL)  in  a 
polypropylene  centrifuge  tube  was  spiked  with  an  approximately  equimolar  amount  of 
f’sZn.  After  several  hours  of  isotopic  equilibration,  samples  were  precipitated  with  300  pL 
of  Optima  trace-element  free  ammonium  hydroxide  (Fisher  Scientific),  The  precipitate 
was  concentrated  by  centrifugation  and  fluid  was  poured  off.  The  remaining  precipitate 
was  dissolved  in  1  mL  of2N  1 1C1,  the  C1‘  concentration  at  which  Zn  has  the  highest 
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Figure  5.2.  Micro-column  design  for 
purification  of  Zn  from  Mg(OH)2 
precipitates,  for  measurement  of  Zn 
concentration  in  seawater  by  isotope 
dilution. 
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Figure  5.3.  Comparison  of  the 
measured  Zn  concentration  for  this 
study  and  by  R.  Wisniewski 
(Wisniewski,  2006)  for  a  vertical 
profile  of  samples  from  the  North 
Pacific.  Error  bars  represent  the  2a 
s.d.  for  triplicate  analysis. 
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affinity  for  anion  exchange  resins  (Kraus  and  Moore,  1953).  Samples  were  purified  to 
remove  Mg  and  other  elements  in  the  precipitate  using  similar  anion  exchange 
chromatography  to  that  used  in  the  preparation  of  samples  for  Zn  isotope  analysis 
(Marechal  et  al.,  1999).  75  pL  of  AG  MP-1  resin  (Bio-Rad)  was  loaded  into  a  column 
made  by  inserting  a  small  frit  into  a  large  pipette  tip.  The  column  exterior  was  wrapped  in 
Teflon  tape  and  the  column  was  inserted  tightly  into  an  acrylic  stand  so  that  filtered  N2 
overpressure  could  be  applied  to  quickly  (-5  mL/min)  push  sample  and  eluent  through 
the  column  (Fig.  1).  Samples  were  loaded  in  2N  HC1,  rinsed  twice  with  0.5  mL  2N  HCi, 
and  eluted  with  250  pL  0.06N  HC1.  The  column  was  rinsed  twice  in  between  each  sample 
with  1-2  mL  0.06N  HC1.  Samples  were  analyzed  on  an  IsoProbc  multi-collector  ICP-MS 
using  an  APEX  Q  inlet  system  (ESI)  with  a  75  pL  min'1  MicroMist  nebulizer  (Glass 
Expansion)  (Elemental  Scientific,  Inc.).  Signals  were  monitored  on  masses  64  and  68 
("  'Zn  and  <,hZn)  to  calculate  the  concentration  of  Zn  in  the  original  seawater  sample  by 
isotope  dilution. 

Blanks  were  interspersed  with  all  samples.  Blanks  were  prepared  and  purified  in 
the  same  manner  as  samples,  except  using  was  1  mL  of  Sargasso  Sea  surface  water 
containing  0.33  nM  Zn, 

5.2.4  Separation  of  Zn  from  seawater  for  isotopic  analysis 

The  isotope  effect  of  Mg(OH)2  co-precipitation  and  the  efficiency  ofZn  recovery 
by  this  method  was  tested  by  adding  small  amounts  of  ammonia  to  precipitate  Zn  slowly 
(Table  1.).  Sargasso  Sea  surface  water  (1  L)  was  acidified  to  pH  2,  spiked  with  20  pg  Zn 
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(20  ppb)  and  allowed  to  equilibrate  overnight.  Several  times,  the  entire  sample  was 
precipitated  with  several  mL  of  ammonia,  the  precipitate  was  separated  by  filtration.  The 
filtrate  was  recovered  and  more  precipitate  was  allowed  to  form  by  waiting  a  few 
minutes,  or  by  adding  additional  ammonia.  The  precipitate  from  each  step  was  dissolved 
separately  for  Zn  isotope  and  concentration  analysis. 

A  similar  precipitation  method  was  used  for  standards  and  natural  samples.  Zn 
was  separated  from  seawater  for  isotopic  analysis  by  coprccipitation  with  Mg(OH)2  and 
filtered  to  remove  the  precipitate.  Precipitate  formed  after  the  addition  of  15  mL  Optima 
ammonium  hydroxide  to  1  L  of  seawater.  The  precipitate  was  allowed  to  settle  for  1-3 
hours  and  the  fluid  was  filtered  through  a  0.4  pm  Nuclcporc  polycarbonate  filter 
membrane  (Whatman)  in  a  PFA  filter  holder  (Savillcx).  Approximately  30  mL  of  2N  MCI 
was  used  to  rinse  the  sample  bottle,  dissolve  the  precipitate,  and  rinse  the  filter  holder. 
The  filtrate  was  rcprccipitated  with  5  mL  Optima  ammonium  hydroxide  and  immediately 
filtered  again.  Another  30  mL  of  2N  HC1  was  used  to  dissolve  the  precipitate  and  rinse 
the  labware.  Both  HC1  fractions  were  combined  for  purification  by  anion  exchange 
chromatography. 

5.2,5  Preparation  of  plankton  tows  for  isotopic  analysis 

Plankton  samples  were  scraped  off  the  filter  with  a  clean  plastic  implement  and 
placed  in  8  mL  quartz  beakers.  They  were  reacted  at  low  heat  overnight  with  2mL  16N 
HNO.i  to  bum  off  some  of  the  organic  material.  When  the  sample  was  dry  quartz  beakers 
were  placed  in  a  shallow  beaker  covered  by  a  watch  glass  and  combusted  overnight  at 
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450°C  in  a  muffle  furnace  to  eliminate  residual  organic  material  {Bergquist,  2004), 
Combusted  samples  were  dissolved  in  2N  HCI. 

5.2.6.  Sample  purification  and  isotopic  analysis 

Samples  were  purified  by  anion  exchange  chromatography  on  AG-MPl  resin 
(Bio-Rad)  based  on  previous  methods  (Marechal  et  al.,  1999).  Samples  were  loaded  onto 
the  column  in  2N  HCI,  the  Cl"  concentration  at  which  Zn  has  the  highest  affinity  for  the 
resin  (Kraus  and  Moore,  1953),  in  order  to  conserve  acid  and  maximize  Zn  retention  on 
the  column.  Plankton  tow  samples  were  purified  once  through  a  column  containing 
1,8  mL  of  resin.  Seawater  samples  were  purified  on  a  larger  column,  and  subsequently 
passed  through  a  smaller  heat-shrink  TFE  (Tef  Cap  Industries)  column  with  400uL  of 
resin  to  further  reduce  the  presence  of  contaminates  (Archer  and  Vance,  2004).  Samples 
were  eluted  with  0.06N  HCI  instead  ofO.OSN  HNO^  to  minimize  the  elution  of 
contaminates  (Chapman  et  al,,  2006).  These  purified  solutions  were  evaporated  to 
dryness  in  PFA  capsules,  and  reacted  at  high  heat  overnight  with  400  pL  HNO.i  and  200 
pL  HF  to  remove  residual  silica  from  the  sample  and  residual  organics  that  may  have 
leached  from  the  column.  Isotope  ratios  were  measured  on  an  IsoProbc  according  to 
standard  methods  (Chapters  1-4). 
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5.3  Results  and  discussion 


5.3.1  Zn  concentration  in  seawater 

Isotope  dilution  and  Mg(OH)2  coprecipitation  followed  by  anion  exchange 
purification  has  successfully  been  used  to  measure  Zn  concentrations  in  seawater.  Process 
blanks  contained  an  average  of  0.22  ng  Zn.  corresponding  to  a  blank  of  0.23  nM  for  a 
l5mL  sample.  Standard  deviations  for  triplicate  processing  and  analysis  of  a  single 
seawater  sample  were  between  0  to  0.2  nM  after  discarding  anomalously  high  data  points 
that  were  assumed  to  result  from  sample  contamination.  Less  than  10%  of  the  replicate 
data  points  were  discarded.  Separate  subsamples  from  the  same  GO-FLO  bottles  were 
analyzed  for  dissolved  Zn  concentration  by  R.  Wisniewski  (Wisniewski,  2006). 
Concentration  measurements  agree  well  with  the  exception  of  the  two  samples  with  Zn 
concentrations  above  10  nM  (Fig.  3).  Both  of  these  samples  were  collected  using  a  GO- 
FLO  bottles  attached  to  a  stainless  steel  hydrowire,  instead  of  the  usual  Kevlar.  The 
measured  Zn  concentrations  arc  higher  than  previously  reported  for  the  deep  North 
Pacific,  indicating  that  the  samples  were  contaminated  (Wisniewski,  2006).  The  reason 
for  the  discrepancy  between  our  data  and  that  of  Wisniewski  is  unknown,  but  may  result 
from  the  entrainment  of  small  metal  particles  from  the  hydrowirc  which  caused  different 
levels  of  contamination  in  the  different  subsamples. 

Anion  exchange  purification  reduced  the  amount  of  Mg  in  the  samples  by  over 
two  orders  of  magnitude.  This  reduces  the  possible  24 Mg'10 A r  interference  on  the  t,4Zn 
peak  during  isotope  dilution  analysis.  On  the  IsoProbc,  the  hcxapolc  collision  cell 
effectively  eliminates  most  polyatomic  interferences  and  so  even  unpurified  isotope- 
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Step 


Figure  5.4.  One  liter  of  seawater  spiked  with  20  ppb  Zn  and  precipitated  in  several  steps. 
Instantaneous  Zn  recovery  and  A6^Zn  (♦)  was  measured  for  each  precipitate  to  calculate 
cumulative  Zn  recovery  ( - )  and  the  cumulative  A^Zn  (□)  of  the  recovered  product. 


Step 

Ammonia 

added 

(mL) 

Time  since 
last  step 
(min) 

Zn 

recovered 

(ng) 

Zn  recovery 
(cumulative) 

AS6Zn  of 
recovered 
Zn 

A“Zn 

(cumulative) 

1 

7.5 

2 

12481 

62% 

-0,22 

-0.22 

2 

0 

15 

2322 

74% 

-0,25 

-0.22 

3 

6 

10 

4787 

98% 

0,23 

-O.ll 

4 

10 

60 

816 

102% 

2.01 

-0.03 

5 

10 

60 

48 

102% 

- 

- 

Table  5. 1 ,  Seawater  doped  with  standard  Zn  was  precipitated  in  several  steps.  For  each 
step,  both  the  time  since  the  previous  step  and  the  amount  of  ammonia  added  was 
recorded.  The  amount  and  Zn  isotope  ratio  of  the  Zn  recovered  at  each  step  was  used  to 
calculate  the  cumulative  Zn  recovery  and  the  A^Zn  of  the  cumulative  Zn  recovery. 
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dilution  concentration  samples  can  be  accurately  measured  in  low-resolution  mode.  On 
instruments  without  a  hexapole,  unpurified  samples  must  be  measured  in  a  high- 
resolution  mode  to  eliminate  polyatomic  interferences  (Wisniewski,  2006).  The  sample 
purification  described  here  may  eliminate  the  need  to  make  concentration  measurements 
in  a  high-resolution  mode.  Decreasing  the  amount  of  Mg  in  the  samples  is  also  beneficial 
for  maintenance  of  the  ICP-MS,  After  analyzing  a  large  number  of  samples  with  Mg 
precipitates,  a  white  dusting  of  Mg  can  often  be  observed  on  the  sampling  cone  and 
skimmer  cone.  This  can  decrease  instrumental  sensitivity  and  cause  interferences  (c.g. 
:4Mg4(,Ar)  that  arc  large  enough  to  interfere  with  isotope  ratio  measurements  on  natural 
samples. 

5.3.2  Dissolved  Zn  isotopes  in  seawater 

5.3.2. 1  Zn  separation  from  seawater  by  large-volume  precipitation  and  filtration 

The  recovery  efficiency  of  the  targe-volume  precipitation  method  was  evaluated 
by  precipitating  1L  of  20  ppb  Zn-dopcd  seawater  in  five  steps  (Table  I,  Fig.  4).  Complete 
recovery  of  Zn  was  achieved  by  the  end  of  the  precipitation  Step  4.  The  simplified 
precipitation  process  used  for  actual  samples  is  roughly  equivalent  to  combining  steps  I  -3 
in  a  single  step,  then  completing  a  second  precipitation  equivalent  to  step  4.  Calculated 
recoveries  for  real  samples  using  this  method  was  always  within  error  of  100%. 

The  isotope  effect  of  precipitation  was  determined  by  measuring  Zn  isotopes  in 
precipitate  recovered  from  several  successive  precipitations  of  a  single  seawater  sample 
(Table  1,  Fig.  4).  The  isotope  effect  of  precipitation  was  calculated  for  the  first 
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Figure  5.5.  A  comparison  of  6MZn  and  SM’Zn  measured  in  precipitates  from  a  20  ppb 
seawater  standard  falls  on  the  expected  mass-dependent  fractionation  line  (-  -  -), 
indicating  the  absence  of  isobaric  interferences, 


precipitation  step  using  the  accumulated  product  equation  for  Rayleigh  distillation 
(Mariotti  et  al,.  1 98 1 ).  We  calculate  an  isotope  effect  of  A6f,Zn  -  -0.35%o.  The  preferential 
incorporation  of  light  isotopes  into  the  precipitate  suggests  a  kinetic  isotope  effect. 

The  quantity  and  isotopic  composition  of  the  Zn  recovered  at  each  step 
was  used  to  calculate  the  cumulative  b'^Zn  of  recovered  Zn.  The  isotopic  composition  of 
total  Zn  recovered  was  A66Zn  =  -0.03%o.  Within  our  levels  of  analytical  uncertainty,  there 
is  no  isotopic  offset  between  the  Zn  originally  added  to  the  seawater  and  the  cumulative 
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666Zn  (%0) 


Figure  5.6.  Exponential  mass-dependent  fractionation  will  produce  a  linear  relationship 
between  66SZn  and  6w‘Zn  with  a  slope  of  1.985  ( — )  and  isobaric  interferences  will  cause 
a  point  to  appear  off  the  mass-dependent  line.  The  direction  in  which  the  measured  value 
will  deviate  from  the  true  value  is  shown  by  the  arrows  for  an  interference  on  mass  64, 
66,  and  68. 
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Zn  recovered  from  the  seawater.  Comparison  of  the  6(,('Zn  and  8<lhZn  values  indicate  the 
absence  of  isobaric  interferences  (Fig.  5) . 

5.3. 2.2  Isotopic  analysis  of  natural  seawater  and  standard  seaw  ater  with  low  Zn 
concentrations 

At  the  Zn  concentrations  found  in  natural  seawater  samples,  analytical  problems 
were  discovered  due  to  isobaric  interferences  on  one  or  more  of  the  Zn  isotopes. 
Concentrated  and  purified  natural  seawater  samples  analyzed  here  had  Zn  concentrations 
from  0.7  to  8.6  nM,  equivalent  to  roughly  0.05  to  0.5  ppb.  Seawater  standards  were  also 
processed  with  0.1,  0.5,  and  1.0  ppb  Zn  added.  While  isobaric  interferences  are  rarely 
observed  in  our  other  studies  ofZn  isotopes,  including  the  seawater  experiments  with  20 
ppb  Zn,  the  low  concentrations  of  Zn  in  our  natural  samples  may  increase  the  possibility 
of  matrix  interference. 

The  presence  of  isobaric  interferences  was  ascertained  by  comparing  values  of 
668Zn  and  b^Zn.  Mass  dependent  fractionation  should  produce  a  linear  relationship 
between  6<>sZn  and  666Zn  values  with  a  slope  of  about  two.  Deviations  from  the  mass 
dependent  fractionation  line  indicate  an  isobaric  interference  with  the  same  mass  as  one 
or  more  of  the  Zn  isotopes  (MZn,  ^’Zn,  and  (,KZn).  In  a  two-dimensional  plot  of  bf,tiZn 
versus  666Zn,  the  anomalous  measurement  will  lie  off  the  mass-dependent  line.  The 
direction  of  the  deviation  is  determined  by  the  masses  on  which  the  interference  occurs. 
For  interferences  on  a  single  mass  (64,  66,  or  68),  the  direction  of  the  deviation  is  shown 
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A 


Zn  added 

SMZn  (%o) 

2  a 

s.d. 

66BZn  (%o) 

2o 

s.d. 

0.1  pg 

-0.45 

0.18 

-0.22 

0.14 

0.1  pg 

-0.40 

0.01 

-0.27 

0.11 

0.5  pg 

-0.09 

0.08 

0.42 

0.55 

1  P9 

-0.01 

0.21 

0.09 

0.42 

1  pg 

0.02 

0.10 

0.17 

0.18 

B 

Location  Depth 

[Zn] 

566Zn  (%»] 

2a 
i  s.d. 

2  a 

fi68Zn(%,)  s.d. 

N.  Pac 

26 

25 

-0.01 

0.03 

0.57 

0.14 

it 

60 

25 

0,17 

0.24 

0.73 

0,34 

h 

100 

25 

-0.16 

0.21 

0.26 

0.07 

*i 

150 

50 

0,06 

0.08 

0.18 

0.33 

H 

150 

25 

0.18 

0.07 

0.60 

0,37 

*1 

300 

50 

0,53 

0.06 

1,47 

0.13 

*1 

300 

25 

0,58 

0.12 

1.19 

0.37 

II 

1000 

50 

0.62 

0.08 

1.45 

0.14 

It 

1000 

50 

0.52 

0.06 

1.35 

0.19 

It 

1000 

25 

0,46 

0.13 

0.69 

0.24 

It 

1000 

50 

0.41 

0.02 

0.92 

0.23 

It 

1000 

25 

0.50 

0.15 

1.26 

0.23 

N.  Atl. 

3000 

50 

0,48 

0.05 

1.22 

0.14 

tl 

3000 

25 

0,45 

0.03 

1.14 

0.35 

Table  5.2.  S^Zn  and  &66Zn  values  measured  for  seawater  doped  with  0. 1 , 0.5,  or  1  ppb  Zn 
(  A)  and  for  natural  seawater  samples  (B).  Errors  arc  reported  for  triplicate  analysis  of  a 
single  sample. 
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Figure  5.7.  6hSZn  and  6^Zn  values  measured  for  seawater  doped  with  0. 1 ,  0.5.  or  1  ppb 
Zn  (A)  and  for  natural  seawater  samples  (B).  The  deviation  from  the  line  of  mass- 
dependence  (-  -  -)  indicates  the  presence  of  isobaric  interferences. 

in  Figure  6. 

We  can  see  from  a  plot  of  &68Zn  versus  &f,f’Zn  measurements  of  five  seawater 
standards  prepared  from  spikes  seawater  that  interferences  arc  occurring  on  more  than 
one  mass  (Fig  7A).  Of  the  three  seawater  standards  for  which  anomalous  isotope  ratios 
were  measured,  two  appear  to  deviate  from  the  expected  value  of  A<,(,Zn  =  0  in  the 
direction  expected  for  an  interference  on  mass  64,  the  third  anomalous  value  is  deviant  in 
the  direction  consistent  with  an  interference  on  mass  68.  Multiple  simultaneous 
interferences  are  also  possible. 
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Mass  interferences  also  appear  in  our  measurements  of  natural  seawater  Zn 
isotopes  (Fig.  7B).  With  these  samples,  too,  we  tested  for  the  presence  of  a  single-mass 
interference.  The  correction  for  a  single-mass  interference  can  be  represented  graphically 
by  drawing  a  line  through  the  observed  data  point  with  a  slope  shown  in  the  arrows  in  Fig 
6.  The  8w,Zn  value  at  which  this  line  intercepts  the  mass-dependent  fractionation  line  is 
corrected  for  the  interference. 

Seawater  samples  from  150  m,  300  m,  and  1000  m  in  the  north  Pacific  and  from 
3000  m  in  the  North  Atlantic  were  analyzed  during  multiple  analytical  sessions.  Because 
the  most  common  isobaric  interferences  occur  at  mass  64,  each  sample  was  corrected  for 
the  deviation  from  the  mass-dependence  line  based  on  the  assumption  that  the 
interference  was  on  mass  64  alone.  Isobaric  interferences  on  mass  68  can  be  accounted 
for  by  comparing  these  values  to  the  86e,Zn  value  which  takes  account  of  only  <l(,Zn  and 
^Zn.  w’Zn  interferences  were  not  examined  because  mass-dependent  anomalies  in  this 
direction  (to  the  right  of  the  mass-dependence  line)  were  not  observed.  Correcting  our 
natural  seawater  6M'Zn  values  for  a  mass  64  interference  in  this  way  docs  not  improve  the 
reproducibility  of  our  data  (Fig  8).  Based  on  this  result,  we  presume  that  there  arc 
multiple  isobaric  interferences  in  these  measurements. 

The  source  of  these  interferences  is  currently  unknown.  Only  Ni  has  stable 
isotopes  on  the  same  masses  as  Zn,  and  wNi  corrections  arc  insignificant  to  the  wZn 
measurements,  so  the  interferences  must  be  from  polyatomic  ions.  These  interferences 
may  result  from  incomplete  purification  of  the  samples.  Seawater  is  a  complicated  matrix 
with  a  total-salt/Zn  ratio  of  1 0V- 1 0 1 1  so  even  a  small  remnant  of  the  original  contaminant 
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Figure  5.8.  666Zn  values  for  seawater  samples  measured  during  different  analytical 
sessions.  Samples  are  from  150  m  (A),  300  m  (B),  and  1000  m  (C)  in  the  North  Pacific, 
from  3000  m  in  the  Atlantic  (D)  and  Zn-free  seawater  doped  with  varying  concentrations 
of  standard  Zn  (E).  Shown  arc  the  6w,Zn  values  measured  with  no  correction  for  isobaric 
interferences  (O)  and  corrected  for  an  interference  on  mass  64  (♦).  This  correction  docs 
not  improve  the  reproducibility  of  the  measurements.  Each  point  (per  depth)  represents  a 
separate  analytical  session  or  separate  sample. 
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Figure  5.9.  Dissolved  6w’Zn  values  at 
several  depths  in  the  north  Pacific  and 
north  Atlantic  oceans  (♦),  Error  bars 
represent  the  minimum  external 
analytical  error  (2s  s.d.).  A  trend  towards 
heavier  bb('Zn  values  in  deeper  water 
correlates  with  an  increase  in  the 
concentration  of  Zn  ( - ). 


Figure  5,10.  A  picture  of  the  processes 
that  could  lead  to  lighter  dissolved  &6flZn 
values  in  the  surface  ocean.  Cells 
preferentially  take  up  light  Zn  to  the  cell 
interior,  while  heavy  Zn  is  preferentially 
adsorbed  onto  the  cell  surface.  Cells 
lysis  in  shallow  waters,  and  most  of  the 
internal  Zn  is  released  into  the  upper 
ocean,  Isotopically  heavy  Zn  is 
transported  to  the  deep  ocean  adsorbed 
to  particulate  matter. 
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Figure  5.1 1.  Comparison  of  seawater  dissolved  Zn  isotopes  from  this  study  (♦)  and 
reported  for  another  station  in  the  North  Pacific  (O)  (Bermin  ct  ah,  2006). 


elements  could  cause  interferences.  Interferences  could  also  be  ablated  off  of  the 
instrument.  I  have  observed  that  cleaning  the  sampling  and  skimmer  cones  decreased  the 
magnitude  of  interferences  on  several  occasions.  MgAr  and  SiAr  are  considered  the 
most  likely  interferences.  34Mg4,,Ar  (79%)  interferes  at  mass  64,  26Mg4t,Ar  (1 1%) 
interferes  at  mass  66,  and  2KSi4(,Ar  (92%)  interferes  at  mass  68.  The  fact  that  single- 
isotope  corrections  do  not  improve  the  quality  of  the  data  indicates  that  there  arc  multiple 
interferences. 

53.2.3  Seawater  dissolved  Zn  isotopes 

A  vertical  trend  is  observed  in  the  profile  of  seawater  Zn  isotope  ratios  (Fig.  9).  In 
the  North  Pacific,  6w’Zn  values  are  higher  in  deeper  waters  and  decrease  in  surface 
waters.  The  Zn  isotope  ratio  at  1 000  m  is  the  best  constrained,  with  four  separate  samples 
having  been  analyzed  during  different  analytical  sessions.  The  average  of  all  analyses  is 
6w’Zn  =  +0.50%«  ±  0.08.  The  two  measurements  presented  for  300  m  and  1 50  m  are  for 
the  same  sample  analyzed  twice  during  different  analytical  sessions.  1  believe  the 
observed  decrease  in  6f,f’Zn  at  150  m  compared  to  deeper  waters  is  real.  Two  separate 
measurements  provided  similar  values  of  6M’Zn  and  in  one  of  them  no  isobaric 
interferences  were  observed.  Unfortunately,  samples  at  26  m  and  30  m  were  measured 
only  once  and  strong  isobaric  interferences  were  found  in  both  samples.  Although  larger 
analytical  problems  were  encountered  in  measuring  &M’Zn  in  seawater  samples  from  the 
upper  water  column,  we  believe  the  trend  towards  lighter  66flZn  in  the  upper  ocean  is  real. 
At  1 50m,  two  separate  measurements  confirm  6,,flZn  of  +0.06  to  +0. 1 8%o  and  on  one  of 
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these  measurements,  no  isbaric  interferences  were  observed.  The  seawater  sample  from 
the  deep  Atlantic  has  a  similar  6w,Zn  to  that  of  the  Pacific  samples. 

The  trend  towards  lighter  seawater  8w’Zn  in  the  upper  ocean  is  unexpected  based 
on  the  preference  of  diatoms  for  the  uptake  of  lighter  Zn  isotopes  observed  in  culture 
experiments  (Chapter  4).  In  order  to  produce  isotopically  light  surface  seawater,  Zn 
sinking  out  of  the  surface  waters  must  be  isotopically  heavy  compared  to  surface  waters 
(Fig.  10). 

Contrary  to  conventional  wisdom,  the  Zn  particle  flux  out  of  the  surface  ocean 
may  not  be  “primary”  biological  Zn  in  growing  phytoplankton.  Perhaps  most  of  this 
primary  biological  Zn  is  released  during  cell  lysis  in  the  upper  ocean,  and  the  downward 
Zn  flux  derives  mostly  from  the  adsorption  ofZn:  to  sinking  particles  (Fig.  10).  Zn 
adsorption  to  cell  exteriors  has  repeatedly  been  shown  to  have  a  positive  isotope  effect 
(Chapter  4,  (Gclabcrt  ct  ah,  2006).  While  experiments  indicate  that  Zn  cell-surface 
adsorption  will  be  relatively  unimportant  compared  to  uptake  for  growing  cells  in  the 
surface  ocean,  this  may  not  be  true  of  sinking  particulate  matter. 

To  illustrate  the  difference  between  Zn  and  N  isotopes  during  particle 
degradation,  we  can  consider  the  case  of  a  single  Zn-containing  protein.  When  this 
protein  is  released  from  the  cell  into  seawater  during  lysis,  even  minimal  degradation  of 
the  protein  may  cause  the  Zn  to  be  released.  As  the  protein  degrades  further,  numerous 
charged  groups  will  be  created  and  exposed,  providing  many  sites  onto  which  free  Zn‘* 
atoms  will  rc-adsorb.  Thus  sinking  protein  fragments  may  contain  most  of  thier  original 
atoms  of  N.  while  Zn  atoms  will  have  exchanged  with  Zn  in  seawater. 
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&6eZn  (%o) 

Figure  5.12.  A  comparison  of  and  &66Zn  measured  in  plankton  tow  samples  from 
several  different  sampling  sites  around  the  world.  d^Zn  values  range  from  0.08%o  to 
+0.5  7%o. 


Figure  5.13.  Worldwide  location  of  plankton  tows,  measured  Sf,(’Zn,  and  Zn:P  of  the 
plankton  (pmol/mol), 
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Our  measurements  of  dissolved  seawater  Zn  isotopes  are  significantly  different 
from  measurements  made  at  another  location  in  the  North  Pacific  (48°39’N,  126°40'W) 
(Bcrmin  ct  al.,  2006).  These  differences  may  be  attributed  to  real  differences  in  seawater 
b66Zn  at  these  locations,  or  to  analytical  problems. 

5.3.3  Zinc  isotopes  in  natural  plankton 

5.3.3. 1  World-wide  plankton  6MZn  values 

Zinc  isotope  ratios  in  plankton  from  around  the  world  range  from  +0.09%o  to 
+0.5 \%o  (Fig.  12,  Fig  13).  The  three  heaviest  6<rf,Zn  values  occur  in  the  samples  that  are 
the  farthest  from  the  shoreline  of  a  major  continent.  Conversely,  many  of  the  l  ightest 
6M’Zn  values  occur  close  to  coastlines,  Measuring  plankton  6w’Zn  values  in  a  broader 
range  of  samples  will  help  to  determine  if  this  is  a  common  feature.  Perhaps  the  isotopic 
similarity  between  open-ocean  plankton  and  deep-water  Zn,  and  the  similarity  between 
coastal  plankton  and  anthropogenic/continental  Zn,  represents  differences  in  the  source  of 
Zn  to  these  plankton. 

5.3.3. 2  The  Peru  Upwelling  Region 

Zinc  isotopes  were  measured  in  phytoplankton  collected  at  nine  stations  in  the 
Peru  Upwelling  Region  (Tahlc  3).  Biomass  in  this  area  is  dominated  by  fast-growing 
diatoms  because  of  the  rapid  supply  of  nutrients  to  the  surface  from  deeper  waters.  Much 
of  the  biomass  here  appeared  to  be  diatoms  based  on  the  large  amount  of  silicate  frustulcs 
left  after  acidification  and  combustion  of  the  samples. 
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Unwashed  cells  Washed  cells 


Station 

6“Zn  (%,)  s.d.  (2s) 

66SZn  (%»)  s.d.  (2s) 

6MZn  (%,)  s.d.  (2s) 

6“Z n  (%<>}  s.d,  (2s) 

10 

002 

0.02 

0.17 

0.16 

0,09 

0.03 

0,35 

0,19 

12 

0.44 

0,05 

1.55 

0.14 

0.96 

0.08 

2,15 

0.42 

24 

-4,26 

0.02 

-8,07 

0.10 

-0.31 

0,05 

-0.11 

0.09 

26 

0,07 

0.02 

0.08 

0,39 

-1.18 

0.05 

-2.42 

0.46 

27 

0.02 

0.02 

0,07 

0,02 

0,23 

0.02 

0.46 

0,09 

27 

-0.14 

0.06 

-0,02 

0.09 

0.14 

0.07 

0.54 

0.02 

29 

0,07 

0.05 

0.15 

0.04 

0.12 

0.04 

0.28 

0.09 

30 

0.51 

0.03 

1.18 

0.06 

0.01 

0.03 

0.04 

0.12 

32 

0,09 

0.02 

0.26 

0,05 

0.16 

0.02 

0.41 

0.11 

33 

0.20 

0.03 

0.44 

0.04 

0.22 

0.02 

0.48 

0.06 

B 

Station 

Latitude 

Longitude 

[Si04]  (nM) 

[P04]  (^M) 

10 

-15.6 

-75.1 

33.02 

1.72 

12 

-16.3 

-75.6 

10.36 

0.72 

24 

-12.3 

-79.3 

1,34 

0.91 

26 

-12.0 

-78.7 

1.59 

0.59 

27 

-11.8 

-78.1 

1.50 

0.62 

29 

-118 

-77.9 

7.59 

0.69 

30 

- 

- 

- 

- 

32 

-11.0 

-78.2 

32.82 

1.89 

33 

-8.2 

-80.3 

3.63 

0.28 

Table  5.3.  (A)  Zn  isotope  ratios  of  washed  and  unwashed  phytoplankton  cells  collected  at 
stations  in  the  Peru  Upwelling  Region.  (B)  Station  locations  and  nutrient  concentrations 
in  surface  waters. 
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Figure  5.14.  A  comparison  of  Sf,!iZn  and  6M’Zn  measured  in  plankton  tow  samples  from 
the  Peru  Upwelling  Region.  When  all  samples  that  deviate  from  the  expected  relationship 
between  6w'Zn  and  86!fZn  by  more  than  the  2a  error  are  removed,  the  remaining  samples 
have  6w’Zn  between  0%o  and  +0,3%o. 


Figure  5.15.  The  b66Zn  of  washed  (O)  and  unwashed  (♦)  plankton  tows  from  the  Peru 
Upwelling  Region.  Only  samples  which  showed  the  expected  mass-dependeni 
fractionation  relationship  between  6l,xZn  and  6hl’Zn  were  included.  There  is  no  clear 
relationship  between  nutrient  drawdown  and  dissolved  Zn  isotope  fractionation  in  the 
Peru  Upwelling  region.  Washed  plankton  tows  are  generally  lighter  than  unwashed 
plankton  tows,  though  the  difference  is  small. 
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'  V  ’*  " 

il  P  ^ 

V  f  t 

*  4 

Unwashed 

plankton 

Washed  plankton 

Krill 

Stn  27 

6EBZn  =  +0.20%o 

6BBZn  =  +0.1 9%„ 

6BBZn  =  +0.1 7%o 

Stn  24 

66BZn  =  +0.22%o 

&BBZn  =  +0.20%o 

6bbZn  =  +0.30%o 

Table  5.4.  The  Zn  isotope  ratio  measured  in  natural  samples  of  unwashed,  and  washed 
phytoplankton,  and  in  krill  picked  out  of  the  sample  from  the  same  tow.  Krill  were 
present  an  only  two  stations  in  the  Peru  Upwelling  Region,  Station  27  (upper  row)  and 
Station  33  (lower  row). 


Several  of  the  samples  showed  extremely  light  values  of  S^Zn  (Fig  14A).  The 
stations  at  which  we  recorded  these  anomalously  light  Zn  isotope  ratios  (6h,’Zn  «  0%o ) 
were  next  to  each  other,  and  were  both  oligotrophic.  Perhaps  these  unusual  numbers 
result  from  analytical  problems,  although  the  values  of  6w’Zn  and  S6* Zn  agree  and  the 
instrumental  mass  bias  during  these  measurements  was  comparable  to  other  samples.  If 
these  extremely  low  S^Zn  values  are  real,  perhaps  they  result  of  contamination  with  a 
highly -fractionated  source  ofZn  such  as  the  ship's  sacrificial  anode.  When  all  samples 
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which  diverged  from  the  mass-dependant  relationship  between  568Zn  and  6 66 Zn  were 
removed,  866Zn  for  the  remaining  samples  ranged  between  0%o  and  0,3%o. 

When  these  samples  are  plotted  against  the  concentration  of  silicate,  no 
relationship  is  observed  between  nutrient  depletion  and  Zn  isotope  fractionation  (Fig.  15). 
Assuming  a  model  where  deep  waters  arc  upwcllcd  with  a  full  nutrient  load  and  nutrients 
arc  slowly  depleted  by  biological  uptake,  even  a  small  isotope  effect  for  uptake  should  be 
magnified  by  Rayleigh  distillation.  From  our  data,  it  appears  that  cither  the  isotope  effect 
for  biological  uptake  in  the  Peru  Upwclling  Region  is  very  small,  or  Rayleigh  distillation 
is  limited.  If  nutrient  concentrations  are  governed  by  differences  in  vertical  upwclling 
rates  at  different  locations,  rather  than  by  horizontal  processes,  plankton  would  reflect  the 
b^Zn  of  upwelling  waters  at  steady  state.  This  might  explain  why  there  is  so  little 
variation  in  plankton  8(,l'Zn  values. 

5.3.3.3  Evidence  for  trophic  fractionation 

At  two  stations,  trophic  fractionation  was  assessed  by  comparing  Zn  isotope  ratios 
in  krill  and  the  diatoms  on  which  they  were  grown  (Table  4).  At  neither  of  these  stations 
were  the  krill  significantly  different  from  the  bulk  plankton  on  which  they  were 
presumably  feeding.  Trophic  fractionation  at  these  sites  appears  to  be  too  small  for  us  to 
measure. 
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Chapter  6 


Conclusions 


I  have  reported  measurable  isotopic  fractionations  in  hydrothermal  systems, 
cultured  marine  phytoplankton,  natural  marine  plankton,  and  seawater.  As  more  data  on 
Zn  isotopes  in  the  marine  environment  becomes  available,  the  potential  for  using  Zn 
isotopes  to  trace  marine  biogeochemical  processes  will  grow.  In  order  to  better  interpret 
this  data,  it  is  worth  considering  how  well  we  understand  the  most  basic  matters  of  Zn 
biogcochcmical  cycling  in  the  oceans.  I  present  here  some  speculations  on  the  origin  of 
the  vertical  concentration  profile  of  Zn  in  the  oceans,  the  relationship  between  the 
nutritional  importance  of  trace  elements  and  “nutrient-like”  concentration  profiles,  a 
reconsideration  of  the  factors  responsible  for  the  seasonal  signal  in  sediment  trap  Sw’Zn 
values  observed  by  Marechal  et.  al.  (2000),  and  some  final  thoughts  on  the  factors 
responsible  for  setting  the  6f,(’Zn  of  seawater. 
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6.1  Origin  of  the  marine  vertical  concentration  profile  of  Zn 

The  classic  cxplaination  for  the  concentration  profile  of  nutrients  in  the  ocean  is 
based  on  the  vertical  transport  of  nutrients  {Rcdfield  et  al.,  1963).  As  high-nutrient  waters 
are  upwelled  to  the  surface,  nutrients  are  depicted  by  biological  uptake  into 
phytoplankton.  Nutrients  are  transported  back  to  deeper  waters  by  the  settling  of 
biological  material  that  remineralizes  to  regenerate  the  dissolved  nutrients.  According  to 
such  a  one-dimensional  model,  the  similarities  in  the  Zn  and  Si  profiles  should  arise  from 
similar  rates  of  particulate  remineralization.  However,  Zn  and  Si  do  not  occur  in  the  same 
phases  within  a  cell.  Chemical  dissolution  of  diatoms,  the  predominant  siliceous  marine 
plankton,  show  that  very  little  of  the  cellular  Zn  is  incorporated  into  the  silicate  frustule 
(Collier  and  Edmond,  1984;  Ellwood  and  Hunter,  2000).  X-ray  spectroscopy  has 
confirmed  this  finding  by  showing  that  Zn  is  contained  mostly  in  the  cell  protoplasm  and 
nucleus,  the  same  parts  of  the  cell  where  P  occurs  (Twining  ct  al.,  2003).  It  is  possible 
that  Zn  is  contained  in  a  more  refractory  organic  phase  than  N  and  P,  so  that  Zn 
fortuitously  dissolves  at  a  similar  rate  to  silicate.  Or,  a  one-dimensional  model  may  not  be 
a  good  representation  of  the  processes  that  give  rise  to  the  marine  Zn  profile. 

An  alternative  cxplaination,  based  on  ideas  passed  to  me  from  John  Edmond 
through  Ed  Boyle,  invokes  two-dimensional  mixing  to  account  for  the  similarity  between 
marine  Zn  and  Si  profiles.  Vertical  features  in  the  centra!  ocean  basins  outcrop  along 
isopycnals  as  horizontal  features  in  the  polar  oceans  (Sarmicnto  ct  al.,  2004)  (Fig.  I). 
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Latitude 


Figure  6.1.  A  pictoral  representation  of  the  proccess  by  which  the  biological  depletion  of 
nutrients  in  high-latitude  surface  waters  can  determine  the  vertical  nutrient  concentration 
profiles  low  latitudes.  Zn  and  Si  are  depleted  more  quickly  than  N  in  surface  waters  of 
the  Southern  Ocean,  so  the  simi iarilty  between  the  shape  of  vertical  Zn  and  Si  profiles  in 
the  North  Pacific  could  have  it’s  origins  in  the  surface  Southern  Ocean. 


Because  mixing  along  isopycnals  occurs  much  more  quickly  than  vertical  mixing, 
horizontal  processes  may  be  important  even  if  they  occur  over  much  larger  spatial  scales. 
Rapid  depletion  of  Si  and  Zn,  compared  to  N  and  P,  in  upwelling  Southern  Ocean  waters 
may  explain  nutrient  profiles  in  the  North  Pacific  (Sarmiento  et  al.,  2004;  Sunda  and 
Huntsman,  1 992),  According  to  this  explanation,  the  similarity  in  the  Zn  and  Si  vertical 
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concentration  profiles  may  result  from  similar  uptake  rates  by  diatoms,  instead  of  a 
similarity  in  particulate  remineralization  rates. 

The  isotopic  consequences  of  Zn  biological  uptake  into  phytoplankton  and  Zn 
remineraiization  arc  likely  to  be  quite  different,  lending  hope  that  Zn  isotopes  may  help 
to  clarify  the  relative  importance  of  these  two  processes.  Remineraiization  is  unlikely  to 
have  a  large  isotopic  signal  because  at  steady  state  the  dissolving  Zn  will  have  the  same 
isotope  ratio  as  the  particulate  Zn.  Biological  uptake  by  diatoms  in  culture  has  a 
demonstrated  isotope  effect  of -0.2  to  -0.8  %o  (Chapter  4)  and  the  isotope  effect  may  be 
magnified  during  Zn  depiction  in  surface  waters  by  Rayleigh  distillation. 

6.2  What  causes  some  trace  elements  to  have  nutrient-type  profiles 

6.2.1  What  is  a  “nutrient-type”  profile 

Zn  is  regarded  to  have  a  “nutrient-type  profile  simply  because  the  shape  of  the 
concntration  profile  is  similar  to  that  of  Si,  a  well-studied  nutrient.  Whether  the  vertical 
concentration  profile  of  Zn  is  controlled  by  vertical  processes  alone  or  horizontal  process 
at  high  latitudes,  it’s  nutrient-type  profile  is  often  implied  to  result  from  the  role  of  Zn  as 
a  cellular  nutrient  and  it’s  depletion  to  very  low  levels  in  oligotrphic  surface  waters.  But, 
there  may  be  other  causes  for  a  nutrient- type  profile. 

6.2.2  Cadmium,  intracellular  nutrient  or  extracellular  toxin? 

The  discovery  that  dissolved  Cd  in  the  oceans  has  a  similar  vertical  concentration 
profile  to  N  and  P  (Boyle,  Nature,  1976)  has  led  to  much  subsequent  speculation  about  a 
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biological  role  for  Cd.  In  fact,  Cd  has  been  shown  to  act  as  a  nutrient  to  Zn-stressed 
diatoms  (Price  and  Morel,  1990),  and  Cd-containing  enzymes  have  been  isolated  from 
two  marine  species  (Lane  and  Morel,  2000;  Lane  et  al.,  2005).  Kinetic  data  showing 
increased  Cd  and  Co  uptake  rates  under  conditions  of  low  Zn  is  similarly  ascribed  to  the 
nutritional  effects  of  these  metals  (Sunda  and  Huntsman,  2000). 

In  fact,  a  nutrient-type  profile  may  just  as  easily  derive  from  the  extracellular 
toxicity  of  Cd.  Cd  ions  appear  to  interact  with  the  same  transporters  used  for  uptake  of 
other  metals  such  as  Zn  and  Mn,  and  high  Cd  concentrations  can  retard  cell  growth  by 
interfering  with  transport  of  these  other  metals  (Sunda  and  Huntsman,  1996;  Sunda  and 
Huntsman,  1998;  Tortell  and  Price,  1996).  In  this  case,  Cd  is  toxic  only  when  it  is  outside 
the  cell!  Perhaps  Cd  uptake  is  simply  a  mechanism  to  “detoxify”  Cd  by  removing  it  from 
the  cells  diffusive  boundary  layer  where  it  interferes  with  other  cellular  machinery. 

It  may  not  be  geochcmically  important  in  the  modem  oceans  whether  Cd  uptake 
is  the  result  ofCd’s  role  as  an  intracellular  nutrient  or  as  an  extracellular  toxin.  It  is, 
however,  a  question  that  should  be  addressed  in  seeking  to  understand  the  evolutionary 
pressures  that  lead  cells  to  acquire  metals. 

6.2,3  Is  Zn  uptake  in  the  ocean  governed  by  it’s  role  as  a  nutrient? 

There  is  no  doubt  that  Zn  can  play  an  important  role  as  a  nutrient.  At  low  Zn 
concentrations,  the  growth  of  phytoplankton  can  become  severly  limited  (Sunda  and 
Huntsman,  1992).  Still,  it  is  not  known  whether  the  uptake  ofZn  in  the  ocean  is  driven  by 
a  physiological  Zn  requirement.  The  Zn:P  ratios  of  cultured  phytoplankton  increase  with 
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increasing  Zn  in  the  culture  medium.  This  may  represent  an  increase  in  the  production  of 
Zn-containing  proteins,  but  at  some  point  Zn  uptake  surely  exceeds  the  biological  need. 
Diatoms  grown  in  high-Zn  media  have  been  observed  to  take  up  Zn  at  mole  ratios  of 
nearly  than  1 :1  compared  to  P  (Sunda  and  Huntsman,  1992).  It  is  not  known  at  what  point 
Zn  uptake  by  marine  phytoplankton  exceeds  physiological  use  of  Zn.  It  is  also  unknown 
whether  uptake  of  Zn  beyond  immediate  physiological  needs  represents  (1)  “luxury 
uptake”  with  no  harmful  effects  of  excess  Zn  uptake  and  potential  benefits  under  future 
conditions  of  Zn  starvation,  (2)  “accidental"  Zn  uptake  which  harms  the  cell  but  which 
the  cell  is  incapapable  of  turning  off,  or  (3)  “de-loxifying”  uptake  where  the  cell  gains  a 
benefit  from  reducing  the  extracellular  concentration  of  Zn  in  the  diffusive  boundary 
layer. 

6.3  What  is  responsible  for  the  seasonal  S^Zn  cycle  observed  sediment  trap  material 

A  data  set  with  great  potential  to  help  us  understand  Zn  isotope  cycling  in  the 
oceans  is  a  year-long  record  of  6w’Zn  in  sediment  trap  material  from  off  the  west  coast  of 
North  Africa  (Marechal  ct  al.,  2000)  (Fig.  2).  Biological  uptake  of  Zn  in  the  surface 
ocean  was  suggested  to  be  the  underlying  cause  of  the  seasonal  cycle  in  6w’Zn  values  in 
this  sediment  trap  material.  Depletion  of  light  isotopes  from  surface  waters  by  biological 
uptake  early  in  the  spring  could  lead  to  an  increase  in  surface  water  6w’Zn  later  in  the 
spring  and  summer.  The  S66Zn  values  in  sediment  trap  material  was  thought  to  mirror 
trends  in  surface  waters. 
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Figure  6.2.  The  Zn  isotope  composition  of  sediment  trap  materials  collected  off  the 
western  coast  of  North  Africa  at  the  EUMEL1  site  (18°30’N  2P06'W).  From  (Marechal 
et  at.,  2000). 


Figure  6.3.  Calculated  zinc  flux  with  organic  carbon  ( — ■)  and  continental  dust  ( — )  over 
the  course  of  the  year.  Total  particulate  matter  flux  is  also  plotted  ( — ).  Data  are  from 
the  EUMELI  project  (Morel,  1996),  with  ZmCorg  estimated  from  plankton  tow  data 
(Collier  and  Edmond,  1984),  and  Zn:Al  estimated  from  measurements  of  African  dust 
(Reid  et  al„  2003). 
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Interpretation  of  this  data,  however,  may  not  be  as  straightforward  as  originally 
proposed.  After  estimating  the  amounts  of  Zn  arriving  in  the  sediment  traps  with  different 
phases,  we  found  that  changes  in  the  source  of  Zn  provide  an  alternative  explanation  of 
the  observed  seasonal  trends.  Sediment  trap  material  recovered  from  the  EUMF.LI 
program  is  a  mix  of  different  phases  (Morel,  1996)  which  may  contain  Zn.  The  relative 
amounts  of  different  phases  changes  with  the  seasons.  During  the  spring  and  summer, 
there  is  a  greatly  increased  flux  of  biological  carbon  into  the  sediment  traps,  w'hile  dust 
flux  (as  measured  by  Al)  is  much  less  variable.  We  estimated  the  amount  of  Zn  delivered 
to  the  sediment  traps  with  these  two  important  Zn-bearing  phases.  The  Zn  delivery  with 
organic  carbon  was  estimated  based  on  a  Zn:Corg  ratio  of  1 .8*  1 0^  (mol/mg)  in  planton 
tows  from  the  Central  Tropical  Pacific  (Collier  and  Edmond,  1984).  The  Zn  flux  with 
continental  dust  (as  measured  by  Al)  was  estimated  based  a  Zn:Al  ratio  of  1.8*10 5  for 
African  dust  (Reid  et  al.,  2003).  According  to  our  calculations,  there  is  a  dramatic 
increase  in  the  flux  of  Zn  with  organic  carbon  during  the  spring  and  summer,  at  the  same 
time  as  there  is  a  marked  increase  in  the  5^  Zn  of  sediment  trap  particles  (Fig  3). 

Based  on  this  new  interpretation  of  the  data,  we  can  infer  values  of6”"Zn  for  dust- 
bome  and  organic-bom  Zn.  The  6w,Zn  of  sediment  trap  material  collected  during  the 
winter  appears  to  be  influenced  mainly  by  continental  dust  alone,  with  S^Zn  -  +0.2%<>. 
This  value  is  remarkably  consistent  with  the  isotope  value  measured  for  Nigerian  dust  of 
6(,l’Zn  =  +0.1 9%o  (Marechat  et  al,,  2000),  Our  calculations  suggest  that  both  dust  and 
organic  carbon  are  responsible  for  transporting  Zn  to  the  sediment  traps  during  spring  and 
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summer.  Based  on  this,  we  suggest  that  the  spring/summer  6w’Zrt  of  between  +0.3%o  and 
+0.4%o  represents  a  minimum  S6flZn  for  the  Zn  flux  associated  with  organic  carbon. 

The  region  where  this  sediment  trap  was  located,  just  to  the  west  of  the  Saharan 
desert,  has  some  of  the  highest  dust  fluxes  in  the  world  (Jickclls  ct  al.,  2005).  Sediment 
trap  samples  from  elsewhere  might  provide  a  better  record  of  the  biological  flux  ofZn 
isotopes  into  the  deep  ocean, 

6.4  What  controls  the  Zn  isotope  ratio  in  the  deep  ocean? 

1  propose  that  the  Zn  isotopic  composition  of  crustal  material  is  b66Zn  =  +0.2%o  ± 
0.  l%o.  This  is  similar  to  the  isotopic  composition  of  basalt  Zn  {see  Chapter  3).  Most 
anthropogenic  Zn  samples  have  bfi6Zn  =  +0.1  %o  to  +0.3%o,  also  similar  to  this  proposed 
crustal  isotope  ratio  {Chapter  2).  The  isotope  composition  of  the  deep  ocean,  however,  is 
about  b^Zn  -  +0.5%o.  There  must  be  a  large  marine  sink  for  light  Zn, 

Based  on  the  limited  data  available,  many  sinks  for  Zn  in  the  ocean  will 
preferentially  sequester  heavier  isotopes.  We  suggest  that  the  flux  of  Zn  out  of  the  surface 
ocean  is  isotopically  heavy  compared  to  surface  waters  (Chapter  5).  Sediment  trap  data 
support  this  assumption  by  showing  that  biogenic  particles  are  at  least  6M’Zn  =  +0.3%o  to 
+0.4%o.  Adsorption  onto  sediments  or  other  solid  phases  is  likely  to  sequester  heavy  Zn. 
Zn  adsorption  onto  cell  cultures  (Chapter  4)(Gelabert  et  al.,  2006)  and  manganese 
nodules  (Marcchal  ct  al.,  2000)  arc  associated  with  a  positive  isotope  effect. 

The  only  process  we  have  discovered  that  might  create  a  sink  for  light  Zn  in  the 
ocean  is  the  hydrothermal  precipitation  of  Zn  sulfides.  Further  research  may  establish 
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Figure  6,4.  Comparison  of  measured  seawater  {♦)  with  the  6w’Zn  predicted  by  our 
model  ( - )  based  on  isotope  mixing  and  measured  Zn  concentrations  ( — ). 
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whether  this  is  the  process  responsible  for  setting  the  6fihZn  of  the  deep  oceans,  or  we 
may  discover  new  processes  that  control  zinc  isotope  concentrations  in  the  oceans. 


6.5.  Modeling  the  seawater  Zn  isotope  profile 
6.5.1.  Isotope  mixing  mode) 

Based  on  the  hypothesis  developed  in  Chapter  5,  that  Zn  transport  to  deeper 
waters  occurs  primary  as  Zn2’  adsorbed  on  sinking  particles,  a  simple  model  can  be 
constructed  to  explain  the  observed  North  Pacific  Zn  isotope  profile  in  terms  of  isotope 
mixing  between  surface  and  deep  waters  plus  dissolution  of  particulate  Zn.  In  this  model, 
we  specify  Zn  fluxes  and  6w'Zn  for  the  surface  ocean  and  the  deep  ocean  in  order  to 
calculate  66ftZn  at  intermediate  Zn  concentrations  as  an  isotope  mixture  of  surface  and 
deep  waters.  We  assume  that  the  Zn  flux  from  deep  waters  to  surface  waters  occurs  by 
advcction-diffusion,  and  that  Zn  transport  from  surface  waters  to  deeper  waters  occurs 
via  adsorption  to  sinking  particles.  The  deep  water  6w’Zn  is  set  at  66(’Zn  =  +0.50  %o,  the 
average  of  five  measurements  at  1000  m.  Because  we  assume  that  Zn  flux  downwards 
from  the  surface  ocean  is  entirely  due  to  particle  flux  (rather  than  mixing),  the  6w’Zn  of 
the  particle  flux  must  also  be  b^Zn  =  +0.50  %o.  At  intermediate  depths  between  the 
surface  and  deep  ocean,  any  Zn  in  excess  of  the  surface  Zn  concentration  is  assumed  to 
result  either  from  mixing  with  deep  water  (S66Zn  =  +0.50  %a)  or  dissolution  of  sinking 
particles  (6b(,Zn  =  +0.50  %o).  Because  both  of  these  have  the  same  isotope  composition, 
the  dissolved  seawater  666Zn  can  be  calculated  as: 


^  seawater  ~  ^  ^surface  "  / surface  +  & 


deep  /  particulate  J  deep  t  particulate 


fa, 
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where  f surface  is  the  fraction  ofZn  at  that  depth  from  surface  water  and  fdeep/panicuime  is  the 
fraction  coining  from  mixing  of  deep  waters  or  particle  dissolution. /is  calculated  using 
Zn  concentration  measurements  from  Chapter  5.  Though  this  model  is  quite  simple  and 
the  data  are  sparse,  isotope  mixing  modeled  values  of  d66Zn  fit  the  observations 
adequately  (Fig.  4).  The  observed  &fi6Zn  for  surface  waters  in  the  North  Pacific  (at  26  m) 
is  -0.01  %o,  so  adsorptive  isotope  effect  for  this  simple  model  would  be: 

=  6<*Znpor,iei*s  ~  =  +  0.5 1  perm'll 


6.5.1.  Box  model 

A  five-box  model  of  the  upper  ocean  was  constructed  in  order  to  explore  in  more 
detail  the  processes  that  may  lead  to  the  Zn  isotope  profile  observed  in  the  North  Pacific. 
This  model  specifies  a  flux  of  dissolved  Zn  between  each  of  the  boxes  dependant  on  the 
water  flux,  a  flux  of  “primary"  biological  Zn  formed  in  the  surface  ocean,  and  a  flux  of 
“secondary"  Zn  that  adsorbs  to  the  surface  of  particles  (Fig.  5  A).  Dissolved  and 
particulate  fluxes  ofZn  between  the  boxes  were  specified  by  the  equations: 


Box)  (surface  mixed  layer): 

1/A  *7 

— =  -  Q- Zn,  +  Q- Zn, -FpKa- Fp'v-K'-Zn,  *  a 
AT 

Box  2: 

=  -  2  Q  -Zn,  +  Q  -Zn,  +  Q  Zny  +  Fp  -  K  a  +  Fp\2K'Znt  ■  a  -  Fp\y  K'Zn ,  •  a 
AT 


Box  3: 

~Z>i-  =  -  2  -QZn3  +  QZn2  +  QZn4  +  Fp\ a  K'  Zn,  ■  a  -  Fp'u  -K'  Zn,  ■  a 
AT 
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Figure  6.5.  (A)  A  schematic  of  the  5-box  model  showing  water  fluxes  carrying  dissolved 
Zn  (Q*Zn),  the  flux  of  primary  biological  Zn  out  of  the  surface  ocean  (Fp),  and  the  fluxes 
of  particulate-adsorbed  secondary  Zn  (Fp1).  (B)  The  results  of  a  run  of  the  box  model 
where  Fp  =  0.  and  Fp’12  draws  Zn  directly  to  the  deep  ocean  (box  5)  instead  of  box  2. 
This  represents  a  water  column  in  which  particles  formed  in  the  surface  ocean  are 
refractory  and  don’t  remineralize  before  they  reach  the  deep  ocean.  (C)  Results  of  a 
model  run  where  Fp  =  0,  and  the  flux  of  Zn  with  Fp’12  goes  is  split  equally  between 
boxes  2  through  5,  instead  of  going  entirely  to  box  2.  This  is  a  more  realistic  treatment  of 
particle  remineralization,  where  particulate  Zn  is  continuously  dissolving  as  it  sinks  from 
the  surface  ocean  to  the  deep  ocean. 
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Figure  6.5.  (continued)  (D)  Results  of  a  model  run  where  Fp  =  0  and  Fp'  is  the  same  for 
all  boxes.  Because  adsorption  of  Zn  to  falling  particles  decreases  with  Zn  concentration, 
Zn  removal  is  greatest  from  deep  waters  and  the  Zn  concentration  profile  develops  an 
unusual  “concave”  shape.  (E)  Results  of  a  model  run  where  Fp  =  0  and  Fp’  decreases  by 
a  factor  of  five  for  each  box  deeper  in  the  ocean.  By  assuming  that  the  particle  flux 
decreases  dramatically  towards  the  deep  ocean,  the  model  is  able  to  adequately  reproduce 
convex  shape  of  both  the  Zn  concentration  profile  and  the  Zn  isotope  profile  that  was 
observed  in  the  North  Pacific,  even  without  the  sinking  and  remineralization  of  primary 
Zn.  (F)  Results  of  a  model  run  where  Fp  and  Fp’  are  both  non-zero,  and  the  isotope  effect 
for  formation  of  primary  biological  Zn  is  negative  (AM’Zn  = -0.5  %o)  while  the  isotope 
effect  for  adsorption  of  secondary  biological  Zn  is  positive  (A6f,Zn  =  +0.5  %o),  This  model 
represents  the  scenario  described  in  Chapter  5,  where  cells  release  most  of  their 
isotopically  light  primary  internal  Zn  in  the  upper  water  column,  and  particulate  Zn  flux 
deeper  in  the  water  column  is  isotopically  heavy  Zn  adsorbed  to  particulate  matter.  This 
model  produces  an  inversion  in  the  upper  water  column  Zn  isotope  profile  similar  to  what 
has  been  reported  in  other  studies  of  Zn  isotopes  in  seawater  (Fig  6). 
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Constant 

Units 

Run  B 

Run  C 

Run  D 

Run  E 

Run  F 

K 

(mol/mol) 

4.2E-06 

6.8E-06 

2.8E-05 

2.8E-05 

1.0E-06 

K' 

(mol/m3)'1 

10 

10 

10 

10 

30 

O 

(m3/yr) 

6.3 

6.3 

6.3 

6.3 

6.3 

V 

(m3) 

50 

50 

50 

50 

50 

a 

K) 

1 

1 

1 

1 

1 

At 

(year) 

0.03 

0,03 

0.005 

0.005 

0.005 

A6eZn 

(%o) 

- 

- 

- 

- 

-0.5 

A66Zn' 

<%») 

0.536 

0.533 

0.14 

0.33 

0.5 

Fp 

(mol  C/yr) 

0 

0 

0 

0 

3 

F’p  12 

(mol  C/yr) 

2.8 

2.8 

3 

3 

3 

F'Pn 

(mol  C/yr) 

0 

0 

3 

0.6 

0.6 

F'P  34 

(mol  C/yr) 

0 

0 

3 

0.12 

0.12 

F’P  45 

(mol  C/yr) 

0 

0 

3 

0.024 

0.024 

Tabic  6.1.  Variables,  units,  and  values  for  five  separate  runs  of  the  box  model. 


Figure  6.  Seawater  dissolved  Zn 
isotopes  measured  in  this  thesis 
and  by  Berm  in  et  al,  (2006),  as 
discussed  in  Chapter  5. 
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Box  4: 
VAZn 

AT 


=  -  2  ■  Q  ■  Zn4  +  Q  ■  Znj,  +  Q •  Zr\b  +  Fp'M  -K'-Zn}  •  a  -  •K'-Zni  •  a 


Box  5  (deep  ocean): 


where  Q  is  the  flux  of  water  between  boxes,  Znx  is  the  concentration  of  Zn  in  box  ,r,  Fp  is 
the  particulate  flux  carrying  primary  biological  Zn,  K  is  the  Rcdficld  ration  ofZn:C  in 
those  particles,  Fp '  is  the  flux  of  particles  carrying  secondary  adsorbed  Zn,  K  ’  is  the  Zn 
adsorptivity  factor  of  particles,  V  is  the  volume  of  the  box  ,  a  is  the  area  on  the  top  and 
bottom  of  the  box,  and  T  is  time.  The  model  was  run  through  one  thousand  interactions 
with  time  step  AT  to  reach  steady-state. 

0  was  derived  as  a  function  of  the  vertical  advective  diffusivity  (/0  of  the  water 
column,  and  the  water  flux  equations  for  each  box: 


(1) 


where  C(,ox  is  the  concentration  of  an  element  in  the  box,  C,op  is  the  concentration  of  the 
element  in  the  box  above,  and  C/,0„0m  is  the  concentration  of  the  element  in  the  box 
below.  Equation  1  can  be  rewritten  as: 


(2) 


Substituting  a'AZ  for  V  and  rearranging  yields: 


AC  Q'  AAC 
AT ~  aAZ 


(3) 
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The  vertical  diffusivity  is  given  by: 


SC  rd2C 
ST  ~  SZ 2 


(4) 


and  the  discrete  rorm  of  this  equation  is: 

AC  AAC 
AT  ~  A Z2 

Combining  equations  3  and  5  yields: 

QAAC  AAC 
aAZ  ~  A Z2 

which  can  be  simplified  to: 


(5) 


(6) 


(7) 


where  AZ  is  the  height  of  the  box  and  a  is  the  square  area  of  the  top  and  bottom. 

An  attempt  was  made  to  choose  oceanographical  ly-relevant  values  for  each  of  the 
model  parameters.  The  Zn  concentration  of  the  deep  ocean  was  set  at  8nM.  and  it’s  Aw'Zn 
was  set  at  +0.5  %o.  A  value  of  10'5  m2  s'1  was  chosen  for  K.  a  typical  value  for  the 
thermocline.  Fp  and  Fp'  were  based  on  estimates  of  POC  export  in  the  North  Pacific 
during  the  summertime  using  Th-234  as  a  tracer  (Charcttc  et  al.,  1999).  The  Zn:C 
Redficld  ratio  was  based  on  plankton  data  from  the  North  Pacific  (Collier  and  F.dmond, 
1984).  In  some  model  runs,  these  parameters  had  to  be  adjusted  to  produce  the  expected 
Zn  concentration  and  Zn  isotope  profiles  (Table  1)  . 

Using  these  values,  the  model  consistently  predicted  a  particle  flux  of  Zn  from 
surface  waters  that  was  too  large  to  be  compensated  for  by  the  vertical  advcctivc- 
diffusivc  flux  ofZn  from  deeper  waters.  Profiles  similar  to  what  we  observed  in  the 
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North  Pacific  could  be  obtained  by  increasing  the  vertical  advective  diffusivity, 
decreasing  the  particle  flux,  or  decreasing  values  of  K  and  K’  (the  amount  ofZn 
transported  with  a  given  particle  flux).  In  order  to  produce  an  isotopic  ofTset  of 
A f ' f 'Z Odcqp-surfacc =  +0.5%o  between  surface  and  deep  waters,  these  models  required 
adsorption  isotope  effects  of  A^’Zn’  =  +0.33  to  +0.53  %o. 

With  this  model,  we  were  able  to  adequately  reproduce  many  of  the  observed 
features  of  the  Zn  concentration  profile  and  Zn  isotope  profile  observed  in  the  North 
Pacific  (Fig.  5B-F),  Experimental  work  on  release  and  rcadsorption  ofZn  and  Zn 
isotopes  during  particle  degradation,  and  further  measurements  of  the  seawater  dissolved 
Zn  isotope  profile,  are  necessary  to  better  constrain  this  model. 
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